ON STABILITY OF SQUARE ROOT DOMAINS FOR 
NON-SELF-ADJOINT OPERATORS UNDER ADDITIVE 
PERTURBATIONS 

FRITZ GESZTESY, STEVE HOFMANN, AND ROGER NICHOLS 

Abstract. Assuming To to be an m-accretive operator in the complex Hilbert 
space H, we use a resolvent method due to Kato to appropriately define the 
additive perturbation T = To + W and prove stability of square root domains, 
that is, 

dom ((T + W) 1/2 ) = dom (T 1/2 ). 

Moreover, assuming in addition that dom (T 1/2 ) = dom ((Tp*) 1 / 2 ) , we prove 
stability of square root domains in the form 

dom((T + Vl/) 1 / 2 ) =dom(T 1/2 ) = dom ((To*) 1 / 2 ) = dom (((T + W)*) 1 / 2 ) , 

which is most suitable for PDE applications. 

We apply this approach to elliptic second-order partial differential operators 
of the form 

-div(aV ■ ) + (§i ■ V ■ ) + div (B 2 ■ ) + V 
in L 2 (f2) on certain open sets CI C R", n G N, with Dirichlet, Neumann, and 
mixed boundary conditions on dCl, under general hypotheses on the (typically, 
nonsmooth, unbounded) coefficients and on dQ. 



1. Introduction 

The principal aim of this paper is to discuss an abstract approach to the problem 
of stability of square root domains of a class of non-sclf-adjoint operators with 
respect to additive perturbations, and to apply it to the concrete case of uniformly 
elliptic second-order partial differential operators of the form 

-div(aV-)+"^ (1-1) 

in L 2 (0) on certain open sets fi C R™, n E N, assuming a E L°°(SY), under various 
boundary conditions on dil, and under fairly general conditions on the (typically, 
nonsmooth, unbounded) coefficients Bj, j = 1, 2, and V. In addition, we discuss in 
great detail the operator 

- div(aV • ) + (JBi • V • ) + div (B 2 ■ ) + V (1.2) 
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in L 2 (R n ), n e N, n ) 2, with optimal T p -conditions on a, Bj, j = 1,2, and V of 
the type 

a e L°°(R ,l ) nx ™, B U B 2 e L n (R n ) n +i°°(R") n , V £ i Tl/2 (R") + i°°(R n ). 

(1.3) 

In a nutshell, we are interested in the solution of the following abstract problem: 
if To is an appropriate non-self-adjoint operator, for which non-sclf-adjoint additive 
perturbations W of To can one conclude that 

dom ((To + W) 1 ' 2 ) = dom (T Q 1/2 ), (1.4) 

and hence conclude stability of square root domains under additive perturbations? 
More precisely, driven by applications to PDEs, we are particularly interested in 
the following variant of this stability problem for square root domains: if To is an 
appropriate non-self-adjoint operator for which it is known that 

dom(T 1/2 )=dom((T*) 1 / 2 ) (1.5) 

is valid, for which non-self-adjoint additive perturbations W of To can one conclude 
that also 

dom((T + W0 1/2 ) =dom(T 1/2 ) = dom ((T *) 1/2 ) = dom (((T + W)*) 1/2 ) . (1.6) 

Of course, in answering such a question one first needs to be able to define the square 
root of a non-self-adjoint operator and it also necessitates a detailed discussion of 
how to interpret the sum of To and W . We postpone a more detailed account of 
these issues for now, and instead describe a bit of the motivation behind posing 
such a problem (see also Remark 2.15 which provides additional motivations in 
terms of trace formulas and Frcdholm determinants). 

First of all, when discussing variational approaches to PDE (eigenvalue) prob- 
lems, formulations in terms of quadratic (resp., sesquilinear) forms have become a 
standard tool and hence the quadratic form of a non-self-adjoint operator, say, T , 
in a complex, separable Hilbert space Ji, naturally enters the scene. Assuming the 
non-self-adjoint operator To to be m-accretive (cf. (2.8)) so its fractional powers 
can be defined in a standard manner, such a sesquilinear form q To , associated with 
To, should be of the type 

q ro (/,5) = ((n) 1/2 fX /2 9)u for f'9 6 dom(q To ), (1.7) 
such that one can expect to have 

1t (/,s) = (f,T g)n for / e dom(q To ) and g e dom(T ). (1.8) 
However, this instantly begs the question of what to consider for dom(q T )? Given 
(1.7), it would be natural to hope that the common square root domain in (1.5) 
can serve for dom(q T ), that is, one would hope for 

dom(q T J = dom (T 1/2 ) = dom ((T*) 1 / 2 ), (1.9) 

and hence it is natural to hope for the same in connection with its analog for the 
perturbed situation with To replaced by To + V, see (1.6). 

As it turned out, for an m-accretive operator To, Kato [65] indeed proved in 1961 
that 

dom(T«) =dom((r T), a S (0, 1/2), (1.10) 
and he also showed that equality in (1.10) is violated in general if a > 1/2, but 
the important case a = 1/2, and hence precisely what is needed in connection with 
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sesquilinear forms, was left open at that point. Within a year, a counter example 
(the operator associated with d/dx on Wq ' 2 ((0, oo)) in T 2 ((0, oo))) to relation (1.5) 
was provided by Lions [71] in 1962 (see also [72] for additional material on such 
counter examples). Moreover, a counter example invoking a sectorial operator (cf. 
the paragraph following (2.8)) was provided by Mcintosh [73] in 1972. In particular, 
as a consequence of additional results by Kato [66] , this implies that also the form 
domain and the square root domain differed from each other in Mcintosh's counter 
example. 

Digressing a bit, we note that while these negative results ruled out the validity 
of (1.6) for m-scctorial operators in general, they were not the end of the story. 
Quite to the contrary, and due to a great extent as a consequence of continuing 
efforts by Mcintosh over the following decades, the emphasis shifted to concrete 
PDE situations for which an affirmative answer to (1.5) was feasible, leading to the 
notion of the Kato square root problem for elliptic differential operators. In this 
context we refer to [29], [40], [69], [74], [75], [77], [78], where second-order elliptic 
partial differential operators of the type 

J2 d 3 a hk {x)d k + J2 bj(x)d J+ J2 ^( Cj (x)-)+d(x), (1.11) 

are discussed in L 2 (fl), fl C R™ open, under mild regularity assumptions on dfl 
and the coefficients dj t ki bj, Cj, and d (or certain symmetry hypotheses on these 
coefficients). The problem continued its great attraction throughout the 1990s as 
is evidenced by the efforts in [3], [6], [17], [20], [21], [22], [35], [49], [61], [90], [91], 
culminating in the treatise by Auscher and Tchamitchian [23] (see also [24]). The 
final breakthrough and the complete solution of Kato's square root problem for 
elliptic operators of the type — div(aV • ) on E™ with T "^-coefficients a^fe, that is, 
without any smoothness hypotheses on aj t k, occurred in 2001 in papers by Auscher, 
Hofmann, Lacey, Lewis, Mcintosh, and Tchamitchian [12], [13], [14], [15], [56] (see 
also [7], [54], [55]). For subsequent developments, including higher-order operators, 
operators on Lipschitz domains ft C M. n , L p -estimates, and mixed boundary condi- 
tions, leading up to recent work in this area, we also refer to [8], [9], [10], [11], [25], 
[26], [27], [31], [57], [60], [79], [81, Ch. 8], [92, Sect. 2.8, Ch. 16]. 

Returning to the principal topic of this paper, we will now assume To to be 
an m-accretive operator from the outset and describe a class of non-self-adjoint 
perturbations W of T such that T + W is closed and densely defined (thereby 
giving precise meaning to the "sum" of T and W) and such that the stability of 
square root domains (1.4) holds. To define the sum T of T and W, we use a 
factorization approach W = B* A (permitting a two-Hilbert space approach, where 
A, B act between possibly different Hilbert spaces H and JC) and follow a device 
due to Kato [67], which introduces the resolvent of T rather than T itself in the 
form 

(T-zIn)- 1 = (To-zlu)- 1 

- (T -zI H )- l B*[l K - AiT-zIu)-^^ ^A(T - zI n )-\ 

In particular, under natural hypotheses on A and B, T defined in terms of (1.12), 
represents a closed and densely defined extension of the operator (or form) sum 
T + B* A (cf. Section 2 and Appendix A for details). Given T, and under the 
assumption that T is m-accretive, we then employ the standard representation of 
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square roots of m-accretive operators, (cf., e.g., [68, Sect. V.3.11]) 

1 f°° 

(T + EIu)- 1 ' 2 = - dXX-^iT+iE + X)^)- 1 , E>0, (1.13) 
n Jo 

and combine it with the resolvent equation (1.12) for T in terms of the resolvent 
of To. Using the analog of (1-13) with T replaced by To, and under appropriate 
intcgrability conditions on 

A" 1 \\A(T + (A + E)I H )- 1/2 || B(W)K) \\(T a + {\ + E)I H y^B*\\ B{KHy (1.14) 

with respect to A near +oo, we reduce (T + EI-u)^ 1/>2 to a multiplicative per- 
turbation of (T + EI-n)" 1 / 2 (cf. (2.45)), which enables us to prove (1.4) (and 
subsequently, (1.6) given (1.5)). 

In this context we emphasize that additive perturbation problems have of course 
been studied in the literature. We refer, for instance, to [23], [74], [75], [77], and 
[91], where applications to second-order uniformly elliptic differential operators, 
typically, with bounded coefficients, have been treated. Moreover, the case of three- 
dimensional Schrodinger type operators with (strongly singular) Rollnik potentials 
was studied in [91] . (In addition, we note that multiplicative perturbations were also 
discussed in [16].) Our principal objective in this paper was to develop an abstract 
perturbation approach to this circle of ideas that's sufficiently general to comprise 
applications to PDE operators of the type (1.1) under very general conditions on 
the nonsmooth coefficients a and V . 

Next, we briefly turn to a description of the content of each section: Section 
2 contains our abstract approach to the stability of square root domains (1.6). 
We describe Kato's additive perturbation approach defining the resolvent of T in 
terms of To and W = B* A. and prove our principal stability results, Theorem 2.4, 
Corollary 2.6, and a variation of it, Theorem 2.14, which permits a straightforward 
application to PDE situations in our subsequent Section 3. In it we consider various 
L 2 -realizations of — div(aV • ) and its additive perturbation, an operator of multi- 
plication by the function V, in the Hilbert space L 2 (f2; d n x), n EN. We explicitly 
treat Dirichlet, Neumann, and mixed boundary conditions on <9f2, under varying 
assumptions on 17 (from just an open subset of M™ in the Dirichlet case to f2 C R n a 
strongly Lipschitz domain in the case of Neumann or mixed boundary conditions), 
sec Theorem 3.12. Particular emphasis is also given to the special case f2 = W l 
(cf. Theorem 3.9). Our conditions on V are of the type V € L P (Q) + L°°{Q) with 
p > n/2. The case of critical (i.e., optimal) L p -conditions is the prime objective 
in Section 4 which is devoted to a proof of the fact that the square root domain 
associated with 

- div(aV ■) + (f?! • V • ) + div (B 2 ■ ) + V, (1.15) 

in L 2 (K n ) under the conditions 

a 6 L°°(R") nxI \ B U B 2 G L n (M. n ) n + L°°(R n y\ V E L n ' 2 {W l ) + L°°{R n ), 

(1.16) 

is given by the standard Sobolev space W 1,2 (M. n ) for n ^ 3 (the case n = 2 is 
more subtle, see Remark 4.7 («£)). We also note that in our applications we focus 
primarily on dimensions n > 2. The case n = 1 with C R an open interval will 
be considered separately in great detail in [41]. The principal result of Appendix 
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A establishes, in short, equality between T defined according to Kato's resolvent 
method (1.12) on one hand, and the form sum of To and B* A on the other. 

Finally, we briefly summarize some of the notation used in this paper: Let 7-L 
be a separable complex Hilbert space, (•, •)% the scalar product in % (linear in 
the second argument), and In the identity operator in %. Next, if T is a linear 
operator mapping (a subspacc of) a Hilbert space into another, then dom(T) and 
ker(T) denote the domain and kernel (i.e., null space) of T. The closure of a 
closable operator S is denoted by S. The spectrum and resolvent set of a closed 
linear operator in a Hilbert space will be abbreviated by cr(-) and /?(•), respectively. 
The form sum of two (appropriate) operators To and W is denoted by To + q W. 

The Banach spaces of bounded and compact linear operators on a separable 
complex Hilbert space T-L are denoted by B(H) and Soo(H), respectively; the cor- 
responding ^-based trace ideals will be denoted by Bp(TL), p > 0. The trace of 
trace class operators in % is denoted by tr-^(-), similarly, det^(-) abbreviates the 
Fredholm determinant. The analogous notation B(Xi,X 2 ), B oc> (X 1 , X 2 ), etc., will 
be used for bounded and compact operators between two Banach spaces X\ and 
X 2 . Moreover, X\ =-> X 2 denotes the continuous embedding of the Banach space 
X\ into the Banach space X 2 ■ 

In addition, /„ denotes the n x n identity matrix in C", and we use the conven- 
tion that < abbreviates ^ C for an appropriate constant C > 0. Moreover, the 
symbol « between two norms indicates equivalent norms. Finally, we abbreviate 
LP(ri;d n x) :=LP(n) and i p (f), C"; d n x) :=LP{Q) n . 

2. Some Abstract Results 

In this section we present our abstract results on the stability of square root 
domains: Given an m-accretive operator To, we determine conditions on an additive 
perturbation W such that the square root domain of To + W coincides with that 
of T . 

We start by with an appropriate definition of the sum of Tq and W following 
Kato [67] (see also the slight extension in [42]): 

Hypothesis 2.1. (i) Suppose that Tq is a densely defined, closed, linear operator 
in % with nonempty resolvent set, 

P(T ) ^ 0, (2.1) 

A : dom(A) — > K,, dom(A) C T-L, is a densely defined, closed, linear operator from 
% to fC, and B : dom(£>) — > JC, dom(i?) C %, is a densely defined, closed, linear 
operator from % to K, such that 

dom(A) D dom(To), dom(S) D dom(T *). (2.2) 

(ii) For some {and hence for all) z S p{Tq), the operator — A(Tq — zlu)^ 1 B* , 
defined on dom(i?*), has a bounded extension in K,, denoted by K(z), 

K(z) = -A(To - zlnY^B* e B(JC). (2.3) 

(Hi) 1 g p(K(zo)) for some z 6 p(T ). 

Suppose Hypothesis 2.1 holds and define 

R(z) = (T - zlu)- 1 - (T -zI n )-iB*[I K - K(z)]~ 1 A(T - zl n y\ 

z G {C 6 p(T ) | 1 € p(K(C))}- 
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Under the assumptions of Hypothesis 2.1, -R(z) given by (2.4) defines a densely 
defined, closed, linear operator T in H (cf. [42], [67]) by 

R(z) = (T - zl u y\ z G {C 6 p(T ) 1 1 G P (K(C))}. (2.5) 

Combining (2.4) and (2.5) yields Kato's resolvent equation (in the slightly more 
general form of [42], in which Hq is no longer assumed to be self-adjoint), 

(T - zluT 1 = (T - zly.)- 1 -{T Q -zI H )- l B*[I K - K(z)]~ 1 A(Tq - zI H )~\ 

ze{(e p(T ) | 1 G p(K(C))}. (2.6) 
The operator T defined by (2.5) is an extension of (To + B*A)\^ om (x )ndom(B'A)i 
T 2 (T + Ti*A)| dom ( To ) ndom ( B . yl ). (2.7) 

We add that the operator sum (To + B*^4)|dom(T )ndom(B*yl) can be problematic 
since it is possible that dom(To) n dom(B*A) = {0} (for a pertinent example, see, 
e.g., [85, Sect. 1.6]; see also [88]). In light of (2.7), T defined by (2.5) should be 
viewed as a generalized sum of To and B* A, that is, the perturbation W of To 
has been factored according to W = B* A. Under the additional assumption that 
dom (T 1/2 ) = dom ((T *) 1/2 ), and the additional hypotheses (2.10) on A and B, we 
will show in Theorem A. 3 that T also extends the form sum of To and B* A. 

We recall that a linear operator D in % is called accretive if the numerical range of 
D (i.e., the set {(f,Df)-n G C | / G dom(£>), \\f\\u = 1}) is a subset of the closed 
right complex half-plane. D is called m-accretive if D is a closed and maximal 
accretive operator (i.e., D has no proper accretive extension). One recalls that an 
equivalent definition of an m-accretive operator D in % reads 

(D + CIn)- 1 eB(H), \\(D + W 1 ^^ < [MCT 1 , MO > 0- (2.8) 

One also recalls that any m-accretive operator is necessarily densely defined. D is 
called quasi-m- accretive if for some a G R, D + al^ is m-accretive. Moreover, D is 
called an m- sectorial operator with a vertex 7 G M and a corresponding semi- angle 
6 G [0, 7r/2), in short, m-sectorial, if D is a maximal accretive, closed (and hence 
densely defined) operator, and the numerical range of D is contained in the sector 

In the following, we let S lt g, 7 G K, G [0, 7r/2), denote the sector with vertex 
located at 7 and semi-angle 9, 

&y,9:={CeC||arg(C-7)|<0}, B G [0,vr/2). (2.9) 

It is not a priori clear that T, defined according to (2.5), is m-accretive if To is. 
For this purpose we now state the following sufficient condition of m-accretivity of 
T: 

Lemma 2.2. In addition to Hypothesis 2.1 assume that 

dom(A) D dom(T 1/2 ), dom(B) D dom ((T *) 1/2 ) , (2.10) 
and suppose that To is m-sectorial. If there exist 7 > 0, 6 G [0, tt/2), and < 8 < 1 

||A(T - ^r 1/2 || B(w K) I) (T - z)- 1/2 B*\\ B{KM) < 1 - zeC\5 7 , e , (2.11) 
t/ien T defined by (2.5) is m-sectorial. 
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Proof. It suffices to show that o~(T) is contained in a sector and that T is quasi-m- 
accretive. One observes that 

f(T) C a(T ) U (C\{C 6 C 1 1 e p(K(())}). (2.12) 

On the other hand, one has 

C\S^ 9 c{(eC | \\K(Oh(K) < 1} c {C e C 1 1 e p(tf(C))}. (2.13) 

Thus, 

C\5 7 , 9 c{CeC|iep(X(C))}, (2.14) 

and taking complements in C, 

C\{C e C 1 1 e P (^(C))} c S 7 , e . (2.15) 

Therefore, (2.12) and (2.15) yield 

<t(T) Ca(T )U5 7 , 9 . (2.16) 

Since T is assumed to be m-sectorial, a (To) is contained in a sector. As a result, 
<r(To) U 5 7 g is contained in some sector, say Sy By (2.16), Sy g' also contains 
a(T). 

Since To is quasi-m-accretive there exists a constant t > such that To + tin is 
m-accretive and for which 

S 7 / 1 fl/+tC{CeC|Re(O>0}. (2.17) 

Without loss of generality, we may assume that t > 7. In this case, p(To + tin) H 
p(T + covers the left half-plane, {( G C | Re(C) < 0}, and one concludes that 

\\{(T + (t + z)I H )- 1 \\ m) 

< 1 1 (To + (t + z)I H )- l \\ B(H) + \\(T + (t + z)I n )-iB*\\ B(ICM) 

x ||(/ K -if(-(t + z)))- 1 || f5(K) ||A(T + (t + z)/«)- 1 || B(WK) 

|| (T + (t + »)J«)- 1 || B(7i) + || (To + (t + ^K«)" 1/2 ||e (w) 



x || (T + (t + z)I H )-^B*\\ B{JCin) \\(I K - K(-(t + ^)))- 1 || B(K) 
x ||^(To + (t + ^)I K )- 1/2 || B(WiiC) 

C[Re(z)] _1 , z e C, Rc(z) > 0, (2.18) 

for some constant C = C(Tq,5) > 0. To obtain (2.18), we have used three estimates; 
the first two follow from m-accretivity of To: 

UT + (t + z)I n )- l \L (v , < 7rT~T i \\(To + (t + z)I n )- 1/2 L im «S 



" Re(z)' l|V v ' nj " [Re(z)]V2' 

z e C, Re(z) > 0. (2.19) 

The third is 

1 1 [lie - K(-(t + z))}- 1 ]]^ ^ <T\ zeC,Re(z)>0, (2.20) 
which follows from 

\\K(-(t + z))\\ B[K) < ||A(T + (t + ^r 1/2 || e(?w 



x ||(T + (t + ^)-V2 J B*|| B(;CiH) 
< 1 - S, z e C, Re(z) > 0, (2.21) 
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taking the Neumann series for [I/c — K(—(t + z))]^ 1 into account. One notes that 
(2.21) holds based on (2.11) since Re(z) > guarantees —(t + z) € C\S 7 ,0. In light 
of (2.18), T + tl-n is m-accretive and hence T is quasi-m-accretive. □ 

Hypothesis 2.3. (i) Suppose Tq is m-accretive. 

(ii) Suppose that A : dom(A) — > K,, dom(A) C is a closed, linear operator from 
H to K, and B : dom(B) — > K, dom(B) C H, is a closed, linear operator from % 
to K, such that 

dom(A) D dom (T Q 1/2 ), dom(B) D dom ((T*) 1/2 ), (2.22) 
and that there exist constants R > and Eq > such that 

dX A- 1 \\A(T + (A + EVh)- 1 ' 2 I 
a ) (2.23) 

x ||(r + (A + ^)/ H )- 1 /2 B *|| B( ^ ) <oo, E^E Q , 

poo 

lim / dAA" 1 || J 4(r + (A + S)^)- 1 / 2 || 

x ||(T + (A + £OJ«)- 1/a Si B(KjW) = 0. 
Our principal stability result of square root domains then reads as follows. 

Theorem 2.4. Assume Hypotheses 2.1 and 2.3. If T defined as in (2.5) is m- 

accretive, then 

dom (T 1 / 2 ) = dom (T 1/2 ) . (2.25) 
Proof. As a result of (2.23), one infers that 



and 



Um ||A(T + (A + i?)/w)- 1/2 || B(WK) ||(To + (A + E)/h) _1 / 2 B*|| b ^ w j = 0, 

E ^ E , (2.26) 

and consequently, 



Jim \\A(T + EI n )-V 2 \\ m ^\\(T + EI H )-^B*\\ B{KH) = 0. (2.27) 
If K(z) is defined as in (2.3), the estimate 



\\K{-E)\\ B{iq < ||A(T +i : ;/„)- 1 / 2 || emx) ||(To + i?/«)-i/2s*|| B( ^ ); jB>0 , 

(2.28) 

together with (2.27), shows that 

lim \\K{-E)\\ B(tC) =Q. (2.29) 

In particular, there exists E\ > such that 



l|A'(-£)lkc) < llAiTo + EInr^^^iTo + EIn)- 1 /^*^^ < 1/2, 

E^E-l. (2.30) 

Next, define 

/>oo 

J(£7) = / dXX-^ 2 (T + EI H ) 1/2 {T Q + (E + \)I n )-^B*[I K - K(-(E + A))]" 1 
Jo 

xAiTo + iE + X)!^)- 1 , E^xaa.yi{E a ,E 1 }. (2.31) 
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The integral in (2.31) is to be viewed as a norm convergent Bochner integral. To 
confirm that the integral is indeed norm convergent, one estimates the norm of the 
integrand of (2.31) as follows, 

HA-^CTo + EI H )V 2 {T + (E + \)I H )-iB* 

x [I K - K(-(E + \))}- 1 A(T Q + (E + A)/^)" 1 !!^) 



= A" 1 / 2 1| (T + EI n ) 1/2 (T a + (E + \)I H )- 1 ' 2 {T Q + {E + X)I n )~^B* 

x [I K - K(-(E + A))]-M(T + (E + X)I n )-^ 2 (T + (E + X)I n )- 1,2 \\ B(H) 

A- 1 / 2 1| (To + EI n y/ 2 (T + (E + \)Iu)- 1/2 \\ B(n) 

x \\[I K - K(-(E + A))]- 1 !!^!!^ + (E + \)Ih)- 1/2 \\ B(hx) 



x || (T + {E + X)I H )- 1 ' 2 B*\\ B{KM) || (T + (E + A)^)" 1 / 2 \\ B(Hy (2.32) 

A > 0, E ^ max{T ,Ti}. 

By the iJ°°-calculus for m-accretive operators (cf., e.g., [50, Sect. 7.1.3], [70, Sect. 
11]), there exists Mi > for which 

\\{T + EIh) 1/2 {T + (E + X)I h )- 1/2 \\ b{h) <Afi, A>0, E ^ max{T , Ei}, 

(2.33) 

and by (2.30), 

\\[I K -K(-(E + X))]- 1 \\ B{K) ^2 1 A>0, E^max{E ,E 1 }. (2.34) 
Using (2.33) and (2.34) the estimate in (2.32) can be continued as follows 
WX- 1 / 2 ^ + EI n ) l / 2 {T« + {E + X)I H )-^B* 

x [I K - K(-(E + A))]- : A(T + (E + A)^)" 1 \\ m) 
< 2ATA- 1 / 2 ||A(T + [E + A)/«)- 1/2 || e(W K) 

x ||(To + (T + A)/ H )- 1 /2 B *|| i3(K ^||(T + (S + A)/ H )- 1 / 2 || e(H) , (2.35) 

A > 0, E ^ max{.Eo,Ti}. 
Again by the T/°°-calculus for m-accretive operators one obtains the estimate 

\\(To + {E + X)I n )-^X m ^S^:^ ^ A A < < ^ )) A>0,T>0. (2.36) 

Together with (2.30) and (2.36), the estimate (2.35) yields 
HA-V^To + EI H )V 2 {T Q + (E + \)I n )-iB* 

x [I K - K(-(E + A))]- 1 ^ + (E + A)^)- 1 !!^ 
< A/ 1 S- 1 /2 A -i/2 X(o fl)(A) + 2ATA- 1 X (fl,oo)(A)||^(To + (E + A)/«)" 1/2 || B(WX) 

x \\{T + (E + X)I H )- 1/2 B*\\ B{Kn) , X > 0, E > maxlTo,^}. (2.37) 

By hypothesis (cf. (2.23)), the majorizing function in (2.37) lies in T 1 ((0, oo); dX). 
Thus, the integral in (2.31) is norm convergent, and by (2.37), 



p oo 

\\J(E)\\ B (h) ^2R 1 ' 2 M 1 E- 1 ' 2 +2Ah / dXX- 1 \\A{T + {X + E)I H )- 1 ' 2 \ 

Jr 



B{H,K) 
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x\\(T + (\ + E)I n )-V 2 B*\\ B(W £>max{£o,£i}- (2.38) 
In particular, (2.38) and (2.24) show that 

Urn \\J{E)\\ B{H) =0, (2.39) 

which implies the existence of an E 2 > for which 

|| J(E)\\ B(n) < 1, E>E 2 . (2.40) 

Returning to square roots of m-accretive operators, one recalls the representation 
(cf., e.g., [68, Sect. V.3.11]) 

(T + EIu)- 1 / 2 = - / d\\-^ 2 (T+(E + X)I n )-\ E>0. (2.41) 
n Jo 

Applying the resolvent equation (2.6) to the integrand in (2.41) yields 

i r°° 

(T + EIu)- 1 ' 2 = - / dXX-^^To + iE + ^In)- 1 -(To + iE + X)^)-^* 

x |/ K -X(-(£; + A))]-M(To + ( J B + A)/ w )- 1 ], (2.42) 

E ^ max{E , o,^i,^2}- 



Multiplication of (T + (E + X)Iu)~ 1 B* in (2.42) from the left by the identity 



I n = (T + EI n )- 1 / 2 (T + EI n ) 1 / 2 , E>max{E ,E 1 ,E 2 }, (2.43) 
and use of 

(To + EIn)- 1 / 2 = - dXX-^ 2 (T a + (E + X)I H y\ E > 0, (2.44) 
t Jo 

implies 

{T + EI H y 1/2 = (T Q + EI H )- 1 I 2 [I H --K- 1 J(E)], E^ma X {E ,E 1 ,E 2 }. (2.45) 

Since In - tt~ 1 J(E) is boundedly invertible for E > E 2 (cf. (2.40)), (2.45) proves 
that 

^{{T + EIn)- 1 ' 2 ) =mn((T Q +EI n )-^ 2 ), E ^ max{£ , E u E 2 }, (2.46) 
and consequently, 

dom ((T + EI n )^ 2 ) =dom((T Q + EI n )^ 2 ), E > max{E , Bi, E 2 }. (2.47) 
Since 

dom((T + S/«) 1 / 2 ) =dom(T 1 / 2 ), 

wo ^>0, (2.48) 
dom ((T + EI H fl 2 ) = dom (T 1/2 ), 

relation (2.47) implies (2.25). □ 

Next, we strengthen Hypothesis 2.3 as follows: 

Hypothesis 2.5. In addition to the assumptions in Hypothesis 2.3 suppose that 

dom(T 1/2 ) = dom((T*) 1 / 2 ). (2.49) 

Corollary 2.6. Assume Hypotheses 2.1 and 2.5. If T defined as in (2.5) is m- 

accretive, then 

dom (T 1 / 2 ) = dom (T 1/2 ) = dom ((T*) 1 / 2 ) = dom ((T*) 1 / 2 ). (2.50) 
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Proof. By symmetry with respect to To and T * in our current hypotheses, the proof 
of Theorem 2.4 also implies 

dom((T*) 1 ' /2 ) = dom((T*) 1 < /2 ). (2.51) 

The proof of (2.51) mirrors that of (2.25), replacing T by T * and T by T* (inter- 
changing the roles of A and B) and employing the adjoint resolvent equation (A.l); 
we omit the details. 

The result (2.50) then follows from (2.49), (2.25), and (2.51). □ 

Corollary 2.7. Assume Hypotheses 2.1 and 2.5, relation (2.22), and thatT defined 
by (2.5) is m-accretive. If there exist constants E\ > 0, E2 > 0, C\ > 0, C2 > 0, 
£ > 0, such that 

WAiTo + EIny^W^ ^^CMlnE)- 1 / 2 ^, E > E u (2.52) 



\\(T + EI n )-V*B*\\ Bilcn) ^C 2 (]nE)-V 2 +*, E>E ^ ( 2 . 53 ) 
then (2.50) holds. 

Proof. In view of Corollary 2.6, it suffices to verify (2.23) and (2.24) for some R > 
and some E > 0. (2.52) and (2.53) yield 



X-^Apo + (A + E)I n )-V 2 || B(WiK) || (T + (A + E)I H )~^B* \\ B{KH) 

A" 1 [ln(A + T;)]- (1+2£) (2.54) 

< A" 1 [ln(A)]" (1+2e) , £ ^max^-Ea}. (2.55) 

Since (2.55) belongs to i 1 ((2, 00); dA), relation (2.23) holds with i? = 2 and E = 
max{E 1 ,E 2 }. Moreover, by (2.54), 



dX A" 1 \\A(T + (A + E)I H )- l/2 \\ B{nx) || (T + (A + ES)I H )-^B*\\ B{Kfll) 

J dX A- 1 [ln(A + T)]- (1+2£) , E > max^.Ta}. (2.56) 

Dominated convergence (applied to the r.h.s. of (2.56)) then implies (2.24) (with 
R = 2). □ 

The same strategy of proof also permits one to discuss fractional powers different 
from 1/2. For completeness we briefly discuss this next (but we only focus on the 
analog of Corollary 2.6 since the analog of Theorem 2.4 is clear). 

Corollary 2.8. Assume Hypotheses 2.1 and let a 6 (0,1). Suppose that Tq is 
m-accretive and 

dom (T Q ) = dom ((T*) Q ) . (2.57) 
In addition, assume that A : dom(j4) — > K,, dom(A) C TL, is a closed, linear operator 
from Ti. to K,, and B : dom(£?) — > K,, dom(_B) C %, is a closed, linear operator from 
H to K. such that 

dom(A) D dom (T Q ) , dom(5) D dom ((Tq) 1- "), 

dom(A) D dom (T 1 " ) , dom(B) D dom ((T *) Q ) , 
and there exist constants R > and E > such that 

dX X~ a \\A(T + (A + E)I H )~ a \\ B[HK) (2.59) 
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x \\{T + (\ + E)I H )^B*\\ B{KH) < OO, 



dX \~ a \\ (T* + (A + E)I H )^A*\\ B(lc n) (2.60) 

x \\B{T* + (A + E)I n r a \\ B{njc) < oo, E > E , 

with both integrals converging to zero as E — > oo. If T defined as in (2.5) is m- 
accretive, then 

dom (T a ) = dom (T a ) = dom ((T*) Q ) = dom ((T*) Q ). (2.61) 

Proof. We briefly sketch the basic necessary amendment in the proof of Corollary 
2.6: One now starts with the integral representation for the a-th power of the 
resolvent of an m-accretive operator To (cf., e.g., [68, Remark V.3.50]), 

(T + EI H y a = !^M r dXX~ a (T + (E + \)I n )-\ E>0, ae(0,l). 

I" JO 

(2.62) 

Mimicking the line of reasoning from (2.26)-(2.25), one then uses the a-th power 
analog of (2.26), namely the observation that 



lh^||^(To+(A+T)/ H )- tt || K(W!K) ||(ro + (A + J B)/ w )«-iB*|| s(ycw) =0, E^E . 

(2.63) 

Carefully following (2.27)-(2.25), and appropriately modifying powers line by line 
as necessary, one arrives at 

dom (T Q ) = dom (T Q ) , (2.64) 

having used the set of containments in the first line of (2.58) together with finitcness 
and convergence to zero (as E ~ > oo) of the integral appearing in (2.59) along the 
way. Finally, repeating the same argument with the roles of A and B interchanged 
and To replaced by T *, one proves the equality for the domains of the a-th power 
of the adjoints, 

dom ((T*) Q ) = dom ((T *) Q ), (2.65) 
and (2.57), (2.64), and (2.65) then yield the desired chain of equalities in (2.61). □ 

One notes that because a and 1 — a coincide when a = 1/2, the four domain 
conditions in (2.58) reduce to two conditions when a = 1/2. Similarly, when a = 
1/2, the two integrability conditions in (2.59), (2.60) reduce to a single integrability 
condition, as in this case, a simple argument with adjoints shows the two integrands 
in (2.59) and (2.60) actually coincide. 

Remark 2.9. The assumption in Theorem 2.4 and Corollaries 2.6, 2.7, and 2.8 that 
T) and T are m-accretive may be replaced by the assumption that To and T are m- 
sectorial. If To and T are m-sectorial, To and T are quasi-m- accretive: there exists a 
t > such that T + tIu and T + tl-u are m-accretive. Since the assumptions (2.23) 
and (2.24) on To remain valid with To replaced by To + tin, the proof of Theorem 
2.4 can be carried out with T (resp., T) replaced by T + tl-u (rcsp., T + tl-u) to 
show that 

dom ({T + tin) 1 ' 2 ) = dom ((T Q + tl n y/ 2 ) = dom (((T + f^)*) 1 / 2 ) , (2.66) 
implying 

dom (T 1 / 2 ) = dom (T n 1/2 ) = dom ((T *) 1/2 ) = dom ((T*) 1 / 2 ). (2.67) 
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Lemma 2.10. Let S and T be m-accretive operators in H and suppose that 

dom (T 1/2 ) C dom (S 1/2 ) . (2.68) 
Then there exists a constant C > such that 

sup 11(5 + E) 1/2 (T + E)~ 1/2 \\ Km , <C. (2.69) 

E>0 "O(i-L) 

Proof. Employing the i/°°-calculus for m-accretive operators (cf. [50, Sect. 7.1.3], 
[70, Sect. 11]) once again, there exists a constant M > such that 



sup 

E>0 



(S*) 1/2 (S* + £r 1/2 || B(W) «S M, (2.70) 



sup \\T 1 '\T + E)~V 2 \\ < M, (2.71) 

sup\\EV 2 (S*+E)-V 2 \\ s^M, (2.72) 

sup ll^/^T + S)- 1 / 2 !! s;m. (2.73) 

By (2.68), there exists M' > (cf., e.g., [68, Remark IV.1.6], [82, Corollary 2.6.11]) 
such that 

\\S 1/2 f\\ n < M'WT 1 ' 2 ^, f e dom (T 1 / 2 ), (2.74) 

and, in particular, 

\\S 1/2 (T + E)- 1/2 u\\ H < M'\\T 1/2 (T + E) 1/2 u\\ n , E > 0, u G H. (2.75) 
As a result of (2.71) and (2.75), 

sup ||S 1/2 (T + £)- 1 / 2 |L_ < MM'. (2.76) 

Using (2.70), (2.72), (2.73), and (2.76), one estimates 
\{f,(S + E)^ 2 (T + E)-^ 2 g) n \ 

= | (/, (S + E)- l l 2 {S + E){T + £)- 1/2 . 9 ) w | 
= | ((S* + E)-^ 2 f, (S + E)(T + E)-^ 2 g) n \ 
^ \((S* +E)- 1 ' 2 f 1 S(T + Er 1 ' 2 g) H \ 

+ \{{S* + E)- 1 l 2 f,E{T + E)- 1 ' 2 g) n \ 
= | {{S*fl 2 {S* + E)-^ 2 f, S^ 2 (T + E)-^ 2 g) H \ 
+ \{E 1 ' 2 {S* +E)- 1 ' 2 f,E l ' 2 {T + E)- 1 ' 2 g) n \ 

<||(S*) 1 /2(S* + £)- 1 /2|| s(?f) ||5 1 /a(T + £0- 1/a || B(7i) 

+ ||^ 2 (^ + j e)- 1 / 2 || b(h) ||^ (t + ^)-v 2 || bw) 

^M 2 M' + M 2 , E>0, f,geH,\\f\\ H = \\g\\ u = l. (2-77) 

The estimate (2.69) now follows from (2.77) (cf. [89], Theorem 4.3(c)), choosing 

C = M 2 M' + M 2 . □ 

Corollary 2.11. Let S be a densely defined, closed, linear operator in %. Suppose 
that Tj, j = 1,2, are m-accretive operators in H, with 

dom(Tl /2 ) = dom(T 2 1/2 ) C dom(S). (2.78) 
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Then there exist constants < c ^ C < oo such that 

cpi^ + EI n )~^ 2 \\ < \\S{T 2 +EIu)- 1,2 \\ B(m 

< c||s'(t 1 + j bj w )- 1 / 2 || b( ^ ) , £>0. 

Proof. By symmetry (in T\ and T 2 ), it suffices to show the existence of a C > 
for which the second inequahty in (2.79) holds. Since dom(T 1 1 ^ 2 ) = dom(T 2 1 ^ 2 ), 
Lemma 2.10 provides for the existence of a constant C such that 

\\(T 1 + EI n ) 1 ' 2 (T 2 + EI n )- 1 '' 2 \\ B(;H) < C, E > 0. (2.80) 

Hence, one estimates 

||5(T 2 + J B/ K )- 1 /2|| ^ ||5(T 1 + J B/^)- 1 /2 



B('H) ^ W'-'^t ' IIB(M) 

R7-„.^-V2|| 

B(H) 

TP. S> I 

B("H)' 



x ||(r 1 + J E;/ w ) 1 / 2 (r 2 + £;/ w )- 1 / 2 | 



<C||5(T 1 +£;/ w )- 1 / 2 || KWV £>0. (2.81) 

□ 



Lemma 2.12. Let 5 and T denote non-negative self- adjoint operators in T-L with 
dom(T 1 / 2 ) C dom(5 1 / 2 ). Then items (i) and (ii) below are equivalent: 
(i) There exist constants M\ > 0, Eq 1, and /3\ > smc/i that 

WS^iT + EIuT^W^^MiE-^, E^E . (2.82) 
(n) There exist constants M 2 > 0, £o > 0, and /? 2 > suc/i £/ia£ 

||S 1/2 /|| w <e||T 1 /2/|| w + M 2 e-^||/|| w , < e < e , / G dom (T 1 ^) . (2.83) 
Proof. Assume that item (i) holds. Applying (2.82), one estimates 
ll^/Hi = ||^ 1/2 (T + EI n )-V\T + El u )^ff u 
||^/ 2 (T + EI n )-Wf Bm ||(T + ^) 1/2 /||^ 
MlE~ 2fil [||T 1/2 /||^ + E\\f\\^ 

= M 2 Lr 2/3l ||T 1/2 /||^ + M^ 1-2 ^ 1 Il/H^, E^E , / G dom(T 1/2 ). (2.84) 

Taking square roots throughout (2.84) and using (a + b) 1 / 2 ^ a 1 / 2 + fe 1 / 2 , a, ^ 0, 
yields 

||5 1/2 /|| w < t M 1 E-^\\T 1 ' 2 f\\ n + M l E 1 ^\\f\\ u , E^E , f G dom (T 1/2 ) . 

(2.85) 

Since -E may be taken arbitrarily large, (2.85) shows that S is infinitcsimally form 
bounded with respect to T. Indeed, setting e = MiE^ 131 , (2.85) may be expressed 
as 

||S 1/2 /||„ < eWT^fWu+M^e^WfWu, 0<e< M 1 E^\ f G dom(T^ 2 ). 

(2.86) 

If < ft < 1/2, then (2.86) implies (2.83), choosing 

e = M 1 E^, M 2 = M^\ 2 = iz|^. ( 2 . 8 7) 

zpi 
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If ft ^ 1/2, then one further estimates (2.86) by 

\\S 1/2 .f\\ H < e\\T^ 2 f\\ n + M^e- X \\f\\n, < e < min{l, M r E^}, (2.88) 

/ e don^T 1 / 2 ), 

so that (2.83) and hence item (ii) holds with the choice for M 2 given in (2.87) and 

e„ =min{l,M lJ B -^}, ft = 1. (2.89) 

Next, assume that item (ii) holds. Choosing / = (T+EI^-^g, g e H, E > 0, 
in (2.83), one estimates using the functional calculus, 

\\S^(T + EI n )-^g\\ n 

< e||T 1 / 2 (T + i?/ H ) 1 / 2 ||^ ) +Af 2£ -*||(T + ^)- 1 / 2 || B(H) ]||.g||«, (2.90) 

< [e + M 2S -^E~ 1 / 2 ] \\g\\ H , g e ft, < e < e , E > 0, 
which shows that 

\\S 1 l 2 (T + EI n )- 1 ' 2 \\ B{n) ^e + Mie-^E- 1 ' 2 , < e < e , E > 0. (2.91) 
Choosing e = E~ a for a fixed a satisfying 

< a(l+ft) < 1/2 (2.92) 

and E > e~ 1/q in (2.91) yields 

||SV2( T + £)-l/2|| E>eQ 1,a . (2.93) 

Condition (2.92) guarantees 

E afo-i/2 ^ ^-^ ^ ^ L ( 2 Q4 ) 

Using (2.94) to continue the estimate in (2.93), one obtains 

^^(T + S)" 1 / 2 !!^ < (l + M 2 )E~ a , S>max{l,eo 1/a }. (2.95) 
Thus, (2.82) follows from (2.95), choosing 

M 1 = (l + M 2 ), E =max{l,£o 1/Q }, ft = a. (2.96) 

□ 

Finally, we turn to another stability result which will be applied to concrete PDE 
situations in Section 4. 

Hypothesis 2.13. (i) Suppose that T\, is an m-accretive operator in H with 

dom (Tl /2 ) = dom ((T*) 1/2 ), (2.97) 

and A : dom(yl) — > IC, dom(^4) C %, is a densely defined, closed, linear operator 
from H to IC, and B : dom(_B) — > IC, dom(_B) C 7i, is a densely defined, closed, 
linear operator from T-L to IC such that 

dom(A) D dom (T^ 2 ), dom(B) D dom (T^ 2 ). (2.98) 

(ii) For some (and hence for all) z G p(T\), the operator — A(T\ — zlu)~ l B* , 
defined on dom(_B*), has a bounded extension in IC, denoted by K(z), 



K(z) = -A(Ti - zI n )- x B* € B(IC). (2.99) 



16 



F. GESZTESY, S. HOFMANN, AND R. NICHOLS 



(Hi) 1 6 p(K(zo)) for some zq € p(Ti). 

(iv) Define the densely defined, closed, linear operator T in % by 



(T - zlu)- 1 = (Ti - zluY 1 - (Ti - zIh^B^Ik - K(z)]~ 1 A(Ti - zI n T\ 

ze{(e p{T{) 1 1 6 p(K(C))}. (2.100) 

(v) Suppose T defined by (2.100) is m-accretive. 

Theorem 2.14. Assume that Tq, T\, and T satisfy Hypotheses 2.5 and 2.13. If 

dom (T 1/2 ) = dom (T^ 2 ), (2.101) 

dom (T 1/2 ) = dom (T n 1/2 ) = dom (T} /2 ) 

V ) V o ) \ 1 ) (2 1Q2) 

= dom ((T*) 1/2 ) = dom {{T*) 1 ' 2 ) = dom ((T*) 1 / 2 ). 
Proof. It suffices to verify (2.23) and (2.24) with Tq replaced by Ti, that is, 

/>00 

/ dA A" 1 \\A{T X + (A + B)^)- 1 ^]! 

./« n ^ K) (2.103) 

x ||(Ti + (A + .E)i"- H )- 1 /2.B*|| 



and 

lim / dAA- 1 ||A(T 1 + (A + i?)^)- 1 / 2 || 

e ^°°Jr "B{h,k) (2.104) 



x ||(T! + (A + ^OJ^-VaB*!! w) = 



as (2.102) then follows from Corollary 2.6 with T\ playing the role of To in Corollary 
2.6. We note that Hypothesis 2.13 (in particular, (2.97) and (2.98)) guarantees that 
Corollary 2.6 applies to T\ once (2.103) and (2.104) are verified. 
To verify (2.103), one observes 



\A{T, + (A + T)/«)- 1/2 || s(wx) ||(T 1 + (A + T)/ w )-i/2 jB *|| e(Kw) 
= \\A(T + (A + E)I n )-^ 2 (T Q + (A + E^u) 1 ' 2 ^ + (A + E)I n )- 1 / 2 
x ||{(T X + (A + E)I n )-V 2 (T + (A + E^u) 1 ' 2 x 
x(T + (A + T)^)- 1 / 2 i?*} cl || e(i<;w) 
= II A(T + (A + E)!^- 1 / 2 ^ + (A + E)^) 1 / 2 ^ + (A + E)I H y x l 2 



B(H,K) 

x || [B(T * + (A + E)I n )-V 2 (TS + (A + £)/«) 1/2 x 

x( T i + (A+^r 1/2 ]l B( ^) 

< || A(T + (A + £)/wr 1/2 || 8(WiC) ||B(T * + (A + TKwr 172 !!^^ 

x ||(T + (A + £)/«) 1 / 2 (T 1 + (A + S)/„)- 1 / 2 || 6(w) (2.105) 
x ||(T* + (A + J B)/„) 1 / 2 (T* + (A + T)/ W )- 1 / 2 || B(W) , A>0, £ > 0, 

where for typographical reasons only, the symbol {. . . } cl abbreviates the operator 
closure (which is usually denoted by just a bar in this paper). In view of (2.49), 
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(2.97), and (2.101), Lemma 2.10 guarantees the existence of a constant M > for 
which 

|| (To + EIk) 1 ' 2 ^ + Elny^W < M 1 / 2 , 

, ( ' , E>0. (2.106) 

||(T* + EI n )^ 2 (T* +EI H y 1 l 2 ]\\ B(H) ^ M 1 / 2 , 

Using 



\\BpZ + EI H )- 1 ^ B ^ = \\(T + EI H )-V*B%^ n) , E >Q 1 (2.107) 
and (2.106), one further estimates (2.105) by 



\A(n + (A + E)I H )- l / 2 \\ B(H K) || (3\ + (A + E)I n )-^B*\\ B{K H) 



< Af||A(To + (A + ^)/„)- 1 /2|| B( ^ ) || (To + s/w) -i/2 jB *|| e( ^ ); (2.108) 

A > 0, £ > 0. 

By hypothesis, T satisfies (2.23) and (2.24). With the estimate (2.108), (2.103) 
(resp., (2.104)) follows from (2.23) (resp., (2.24)). □ 

Remark 2.15. We note that the stability of square root domains has important 
consequences for trace formulas involving resolvent differences and symmetrized 
Fredholm (perturbation) determinants as discussed in detail in [44]: Let A and 
Aq be densely defined, closed, linear operators in H- Suppose there exists to G R 
such that A + t I-H and A + t I-x are of positive-type and (A + tolu) £ Sect(wo), 
(Aq + tolu) 6 Sect(o;o) for some ojq G [0,tt). In addition, assume that for some 
ti^ to, 

dom ({A + tJu) 1 ' 2 ) C dom ((A + tilu) 1 ' 2 ) , (2.109) 
dom ((A* +t 1 I n ) 1/2 ) Cdom((A* +t 1 I n ) 1 / 2 ), (2.110) 
(A + hln) 1 ' 2 [(A + hln)- 1 - (A + hl^} (A + hln) 1 ' 2 G Bi(H). (2.111) 
Then 



-ln(det H [{A - Z I n) V* { A - Z I H )- H A - zl n W n (2 n2) 

= tr H ((A - zlu)- 1 - (A - zlu)- 1 ), 



for all z G C\(-i + SU) such that (A - zI- H ) 1 / 2 (A - zI H )~ l {A - zI-h) 1 / 2 is 
boundedly invertible. (Analogous trace formulas can be derived for modified Fred- 
holm determinants replacing Bi(H) by Bk("H), k G N, fc 2, including additional 
terms under the trace.) 

Here A is said to be of positive-type if 

(-oo,0]cp(A) and Ma = sup II (1 + t)(A + tiny 1 II . . < oo, (2.113) 



and A is called sectorial of angle lo G [0,7r), denoted by A G Sect(w), if 

(2.114) 



cr(A) C S u and for all oj' G (w,7r), M(A, u/) = sup \\z(A — zlu) 1 1 1 er-wt < °°> 

zGC\S~7 
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where S u C C, ui G [0,7r), denotes the open sector 

f{zeC\z?0, |arg(z)| <lj}, w6(0,tt), 
[(0,oo), uj = 0, 

with vertex at z = along the positive real axis and opening angle 2u. Assumptions 
such (2.109), (2.110) make it plain that control over square root domains is at the 
core for the validity of the trace formula (2.112). 

In the special case where in addition A and Aq are self- adjoint in %, trace 
relations of the type (2.112) are intimately related to the notion of the spectral 
shift function associated with the pair (A,Aq), and hence underscore its direct 
relevance to spectral and scattering theory for this pair (cf. [44] for details). 

Remark 2.16. We conclude this section by recalling that if q T (-, •) denotes the 
sectorial form uniquely associated with m-sectorial operator T (cf. [68, Theorem 
VI.2.7]), then, as shown by Kato [66], the fact dom (T 1 / 2 ) = dom ((T*) 1 / 2 ) also 
implies 

dom(q T ) = dom (T 1/2 ) = dom ((T*) 1/2 ) = dom(q T „)- (2.116) 

In fact, if two out of dom(q T ), dom (T 1 / 2 ), dom ((T*) 1 / 2 ) are equal, all three are 
equal, see [66], [71], and [81, Theorem 8.2]. 

3. Applications to Divergence Form Elliptic 
Partial Differential Operators 

In this section we apply the abstract formalism of Section 2 to uniformly elliptic 
second-order partial differential operators of the form — div(aV ■ ) + V on certain 
open sets f2 C M n , n <G N, with appropriate boundary conditions on dfi. 

For simplicity of notation, we agree in the following to abbreviate L p (fl; d n x) 
simply by L P (Q), p <E [1, oo) U {oo}, n 6 N. 

We start with some general considerations and eventually strengthen our hy- 
potheses on f2. 

Hypothesis 3.1. Let neN and assume that f2 C M. n is open 1 . 

(i) Suppose a : £1 — > <C nxn is a Lebesgue measurable, matrix-valued function which 
is bounded and uniformly elliptic, that is, there exist constants < a\ ^ 02 < 00 
such that 

allelic- < Re^XaOOo] and |(C,a(x)0c-| < ^IMICIIo, 

f,C G C™ ; a.e. xen. 1 ' ' 

(ii) Denote by W(fi) a closed subspace of W 1,2 (fl) satisfying 

C M (1])CW(!1). (3.2) 

(iii) Let L a ^(n) denote the m-accretive (in fact, m-sectorial) operator uniquely 
associated with the densely defined, accretive, and closed sesquilinear form 



n 



qa,w<n)(/.s) = L* 1 * ((V/)(x),a(x)(V.9)(x)) c , l , f,g £ dom(q aiW(n) ) = W(Q). 

(3.3) 



-^If n = 1, we assume that ft C K equals an open interval. 
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(iv) Suppose that V : £1 — > C is (Lebesgue) measurable and factored according to 

V(x) = u(x)v(x), v(x) = \V{x)\ 1/2 , u(x) = e i& ^ v{x)) v{x) for a.e. x G fl, 

(3.4) 

such that 

W(O) C dom(w). (3.5) 

The quadratic form q^ in L 2 (f2) uniquely associated with V is defined by 

qy(/,g) = {vf,e^^vg) L . 2{ny f,g G dom(q y ) - dom(v). (3.6) 

Under appropriate assumptions on V (see Hypotheses 3.2 and 3.3 below), the form 
sum of q a yy(- ) and q v will define a sectorial form on W(f2) and the uniquely 
associated operator with q a w ^ + q v will be denoted by £ a ,w(fi) +q V ( see also 
the paragraph following (A. 42)). 

With respect to the potential coefficient V we introduce the following hypotheses: 

Hypothesis 3.2. Let n G N ; assume that fi C R™ is open 2 , and let V G L p {Vl) + 
L°°(n) for some p > n/2, p ^ 1 for n ^ 2, and V e L 1 ^) + Z,°°(Q) for n = 1. 

In addition, we also discuss the critical L p -index p = n/2 for V for n ^ 3: 

Hypothesis 3.3. Le£ n G N\{2}, assume that fi C l n is open 2 , and let V G 
L™/ 2 (0) + L°°(fi) /or n > 3 and y G L^fi) + L°°(0) /or n = 1. 

Here U G L q {VL) + L°°(f2) means that V permits a decomposition V — V q + Voo 
with V q G L q (Q) for some q ^ 1 and G 

Remark 3.4. (i) One notes that by assumption (3.1), the real part of q Q w ^ is 
non-negative, 

M*a,vv(n)tfJ)]>ai\\Vf\\Ua), / e W(fi). (3.7) 

(m) To see that the (sesquilinear) form q a w ^ defined by (3.3) is closed, one recalls 
that by definition, q a W r<m is closed if and only if dom(q Q is closed with respect 

to the norm || • \\q a Win) = (Rc[q a W(f2) (-, ■)] + || • ll|2 (0) ) 1/2 - Under condition (3.1), 
II ' Hw^fn) and j| • || q w n are equivalent norms on W(f2). Since W(f2) is assumed 
to be closed under the || ■ ||w 1 . 2 (n) norm, W(Q) is closed under the || • || q w norm 
and it follows that q a w ^ is a closed form. For additional details in connection 
with £ a ,w(f2) associated to (3.3) (especially, in the case where a is a symmetric 
nxn matrix), see, for instance, [5], [45], [46], [47], [48], [51], [52], [53], [59], [62], 
[63], [64], [80], and [81, Sect. 4.1]. 

Remark 3.5. Next, we briefly discuss three special cases corresponding to specific 
choices for W(O) and use this to introduce some convenient notation: 
(i) Choosing W(f2) = 1U ' (il) formally corresponds to imposing a Dirichlet bound- 
ary condition on d£l. For simplicity of notation, wc denote the corresponding op- 
erator L a W i,2^ simply by L a ,f2,D, that is, 

L a ,a,D ■= L a W i,2 {n) . (3.8) 



If n = 1, wc assume that £7 C M equals an open interval. 
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(m) Choosing W(f2) = W ' (fl) formally corresponds to imposing a Neumann 
boundary condition on 90. The corresponding operator L a ,w 1 - 2 (fi) will be denoted 
by L a< a jN , that is, 

L a ,n.N ■= L a ^ w i.2(£i)- (3.9) 
(m) With II denoting a fixed, relatively open subset of dft, consider 



W^(Q) = {/|n | / € C °°(1R™); supp(/) n (dQ\U) = 0} 



■w 1 "'(n) 



(3.10) 



c °°(iR n \(ao\n)) 



where V denotes the closure of T> C W 1,2 (f2) in W 1,2 (f2). (One observes in 

connection with (3.10), that since 9f2\II is relatively closed in dCl, it is also closed 
in R n .) The corresponding operator L a will be denoted by L a ,f2,iii that is, 

L a ,n,n ■= L a ,w£- 2 (n)> ( 3 - u ) 



The choice of W(f2) = Wjj 1 (Cl) formally corresponds to mixed boundary conditions 
with a Dirichlet boundary condition on <90\II and a Neumann-type boundary con- 
dition on n (cf., e.g., [5], [59, Definition 2.3 and Remark 2.4], [64, Definition 2.2], 
[80], and [81, §4.1] as well as the references and discussions therein.). The choice 
Id = is permitted and corresponds to the Dirichlet case (i), 

W^ 2 (Q) = W l 2 {n). (3.12) 

The choice II = dfl is permitted as well, but we caution that in general, Wqq (fi) C 
W 1,2 (ft). However, under additional assumptions, for instance, if ft possesses the 
segment property (cf., e.g., [1, Theorem 3.22]), equivalently, if dfl is of class C (cf. 
[34, Theorems V.4.4 and V.4.7]), it is known that 

Wan (n) - W l >\fi). (3.13) 

Moreover, in the particular case where ft = R™, one infers that W(IR n ) = 
W^ 2 (R n ) by (3.2) and the fact that Cg°(R n ) is dense in W h2 (R n ) (cf., e.g., [1, 
Theorem 3.22]). We then drop the subscript W 1 ' 2 (W 1 ) and simply write 

La '■= L atW l,2( Rn y (3-14) 

Finally, in the special case where the matrix-valued function a coincides with the 
n x n identity matrix /„ a.e. on ft, we shall write 

~ A :=£/„, -Aw(fi) :=£/ n ,w(fi), 

A A j a r (3.15) 

where A := Y^j=i &xj denotes the usual Laplacian on W\ In particular, — 

is self-adjoint, and consequently, its square root domain coincides with its form 

domain, that is, 

dom ((-Awp,) 1 / 2 ) = VV(fi). (3.16) 

Our next result uses the relations (3.17) as input and then proves that if V is 
infinitesimally form bounded with respect to — Ayym), then V is also infinitesimally 
form bounded with respect to L 0j w(n) an d that the square root domains (3.17) 
remain stable under perturbations V satisfying Hypotheses 3.1 (iv). 



STABILITY OF SQUARE ROOT DOMAINS 



21 



(3.19) 



(3.20) 



Theorem 3.6. Assume Hypotheses 3.1 and suppose that £ a .w(n) satisfies 

dom (L$,(n)) = dom ((L* amn) ) 1/2 ) = W(Q). (3.17) 
If there exist constants M > 0, f3 > 0, and Sq > such that for all < e < Sq, 

l|l^l 1/2 /fc (n) ^£||(-A W (n)) 1/2 /||" 2(a) +Me-^||/||i 2(n) , /eW(fl), (3.18) 
fften i/ie following items (i) and (ii) hold: 

(i) V is infinitesimally form bounded with respect to i a ,w(f2) an d 

II l^| 1/2 /||' 2(n) ^Re[q 0iW(a) (/,/)] + Ma7 V^||/||| 2(u) , 

/ G W(fi), < e < a^o- 

(ii) The form sum L aW (n) +q V is an m-sectorial operator satisfying 

dom {(L aiW{Q) + q Vfl 2 ) = dom (((i , w(n) +„ V)*) 1 ' 2 ) 

= dom Kwcn)) = dom (( l :,w W ) 1/2 ) = W(O). 

Proof. Let f)y denote the sesquilinear form associated to the operator of multipli- 
cation by V, that is, 

M/.ff) = Ke iarg(V(0) «ff) L 2 ( n). /,ffGdom(w)=dom(ftv). (3-21) 

Since W(f2) C dom(u) by hypothesis, in order to prove that V is infinitesimally 
form bounded with respect to £ a .w(n)i it suffices to prove the existence of an eq > 
and an 77 : (0, £0) — > (0, 00) for which 

Hll^(n) < ^K W ^)(fJ)} + vmf\\h(n V f e W(fi), < £ < e . (3.22) 

Sufficiency of the condition in (3.22) is due to the obvious estimate 

\W(fJ)\^\\vf\\h (a) , /edomfV), (3.23) 

and the fact f)y is infinitesimally form bounded with respect to q a w ^ if and only 
if the former is infinitesimally form bounded with respect to Re[q a yy(n)] (cf-, e.g., 
[68, p. 319]). 

Freely using (3.18), one estimates 

\\vf\\l Hn} < e\\(-A mQ) )V 2 f\\l 2(n) + Me^||/||| 2(fi) (3.24) 

= e||V/||| w +Me^||/||| 2(n) (3.25) 

<or 1 dte[q 0)W(n) (/ J /)]+Mc-' , ||/||i a(n)j /eW(fl), < e < e 0j 

(3.26) 

where the 2nd representation theorem (cf., e.g., [68, Theorem VI. 2. 23]) was applied 
to the (non-negative, self-adjoint) reference operator — and its associated 

form qj in (3.25). To arrive at (3.26) the first inequality in (3.1) was used. Choos- 
ing e = ai£, < £ < af 1 £o in (3.24)-(3.26) yields the form bound (3.19). 

On the basis of (3.19), £ 0) yy(n) +q V defines an m-sectorial operator by standard 
perturbation theoretic results (see, e.g., [68, Theorem VI. 1.33]). Next, one notes 
that by Theorem A. 3, the operator T defined as in (2.4) and (2.5) in terms of the 
identifications 

7b =ia,w(£J)j A = v, B*=u, (3.27) 
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coincides with the form sum £ a .w(f2) +q V. In particular, according to the factor- 
ization (3.4), L a ,w(f2) an d i a ,w(S2) +q V satisfy the resolvent equation 

(L a ,w(n) +q V - zli,2(p.))~ = (L 0j w(n) - 2^l 2 (o)) _1 



- (io,w(n) - ^i 2 (n)) 1 u[^L2(f2) - ^(^a,w(n) - zI L 2 {n) ) 1 u\ (3.28) 

x "(■C'o,w(n) ~ z-Tr, 2 («)) _ \ « G p( L a,w(n)) n p(£ a ,w(n) +q V"). 

Given (3.27), one now proceeds to verifying items (i), (ni) of Hypothesis 2.1: For 
this purpose one notes that by Lemma 2.12, (3.18) implies the existence of constants 
Af > 0, Ei > 0, and /? > such that 

|K-A W (n)+EJ^ ( n)r 1/2 || B(ia(11)) <M'E-P, E > E 1 . (3.29) 

Additionally, one notes that (3.16), (3.17), and Corollary 2.11 imply that (3.29) is 
equivalent to 

\\v(-A wiQ) +EI L 2 iQ) )- 1 / 2 \\ B{L2{a)) ^M"E-P, E>Ei, (3.30) 

for some constant M" > 0. Subsequently, to verify Hypothesis 2.1 (Hi), one uses 
(3.30) and estimates, 



|| (L a , W (Q) + £7/ i2 (n))" 1/2 «|| B(I ,2 (n)) 
= IK-^wcn) + EI LHn))~ 1/2 \\ B(L 2 m 

= \\ v ( L *a,W(n) + EI L^(n)) 1/2 || B (L2( )) 

< ||w(-A w(n) + EI L 2 (n) )- 1/2 \\ B{L2{n)) 

x ||«(-A W (n) + EI L 2 {n) ) 1/2 (L* a W{n) + EI L 2 (n) )~ 1/2 \\ B(L2[n)) 

< CM'E- , E > E u (3.31) 

where C > is an appropriate constant (guaranteed to exist by (3.16), (3.28), and 
Lemma 2.10) for which 

||«(-A w(n) +EI L 2 [n) f/ 2 {L* a <w{n) +EI L 2 {n) )- 1 / 2 \\ B{L2m < C, E>0. (3.32) 

We note that (3.17) is used freely throughout (3.31). Consequently, (3.30) and 
(3.31) can be used to estimate 



\v{L ayV (n) + EI L 2 {a) ) 1 u\\ B{L2[n)) 



= \\ v {L a ,w(n) + EI L 2^) 1 ^ 2 (L Q ,w(o) + EI L 2( n )) ^Hlg^ 2 ^)) 
^CM'M"E- 2fs , E>E 1 . (3.33) 



Thus, 1 S P(v{L a ,w(n) + EI L 2iq\)- x u) for all E > E\ sufficiently large, implying 
the validity of Hypothesis 2.1 (Hi). 

Thus, by Corollary 2.6 and Remark 2.9, the proof of (3.20) reduces to verifying 
Hypothesis (2.5), that is, Hypothesis 2.3 and (2.49) with the identifications (3.27). 
Hypothesis 2.3 (i) and (2.49) are satisfied by assumption (cf. (3.17)). Regarding 
Hypotheses 2.3(H), (2.22) is clear in light of (3.17) and (3.5). In order to verify 
(2.23), one notes that by (3.30) and (3.31), 



\- 1 \\v(L aiWiQ) + (E + X)I LHQ) )- 1 / 2 \ 



B(L 2 (fi)) 
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B(i 2 (0)) 



x ||(£ a ,w(n) + (E + X)I L 2 {n) )- 1 / 2 u 



< \- 1 CM'M"(E + X)- 213 

< CM'-frfA" 1 " 2 ' 3 , E>E U A > 0, (3.34) 

implying (2.23) with J? = 1, Eq = E\, and the identifications (3.27). An application 
of Lebesgue's dominated convergence theorem, using (3.34), then yields (2.24). 
Finally, the first two equalities in (3.20) follow from Corollary 2.6 and Remark 2.9, 
and the last equality in (3.20) follows from assumption (3.17). □ 

Next we investigate the special case f2 = R™, n £ N, in detail and show how to 
derive the basic inputs (3.17) and (3.18) for Theorem 3.6 in this particular situation. 
We start with (3.18). 

Lemma 3.7. Suppose V satisfies Hypothesis 3.2 with f2 = W 1 , «eN, and let — A, 
dom(— A) = W 2,2 (WL n ), denote the self-adjoint realization of (minus) the Laplacian 
in L 2 (M. n ). Then V is infinitesimally form bounded with respect to —IS. and there 
exists a constant M n ^ p y > (depending only upon n, p, and V) such that for all 
e > 0, 

/ e W 1,a (R n ). 

Proof. Without loss of generality we put the L°°-part Voo of V equal to zero. By 
Hypothesis 3.2 one infers 

«(-A + ^(r))" 1/2 e02 P (£ 2 (l")), E>0. (3.36) 

Since v £ L 2p (R n ), [83, Theorem XI.20] (cf. also [86, Theorem 4.1]) yields 

\\v(-A + EI L2{Rn) )-V 2 \\ B2p{L2iMn)) 

< (2 7 r)-"/^||T/||^ 2 (R „ ) ||(| • | 2 +£)- 1/2 |U (Rn) , E > 0. (3.37) 

Employing 

||(|-| 2 + ^)- 1/a || i ^ (an) =C r n B(' l -*V^, E >0, (3.38) 
where C n > only depends on n £ N, (3.37) and (3.38) imply 

\\v(-A + EI L , (Rn) )- l / 2 \\ B2AL2{Mn)) ^C n , p , v E^- 2 ^l^\ E>0, (3.39) 
where 

C n>p , v := (2^r n ^\\V\\ 1 ^ n) C n . (3.40) 

The form bound in (3.35) now follows by repeating the estimate in (2.84) with 
S = V and T = -A. □ 

Remark 3.8. The infinitesimal form bound result of Lemma 3.7 may be found in 
[36, Proposition 6.4], but the exact form of the bound (in particular, the coefficient 
Af™,p,ye"" /(2p " n) in (3.35)) is not explicitly stated there. The proof given above is 
based directly on this circle of ideas. An alternative proof for the case n 3, in the 
more general setting of domains in R™, may be found in [33, Theorem 1.8.4] (our 
proof of item (Hi) in Theorem 3.14 mirrors this alternative approach). 
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Next we focus on (3.17) in the particular case f2 = R™, n G N. In this context 
we recall that the "Kato square root problem" for second-order divergence form 
elliptic operators in R n , has been solved in [13] (we note that the case n — 1 was 
treated earlier in [29], the case n = 2 in [58]). To be precise, we assume Hypothesis 
3.1(i)-(m) with Q, = R' 1 and hence have W(Q) = W 1,2 (R n ). Then L a defined 
by (3.14) (cf. also (3.3)) formally takes the form of a second-order divergence form 
operator, 

L a := — div(aV • ) = — dja hk dk- (3-41) 

The main result of [13] is that Hypothesis 2.3 (i) holds with T := L , % = L 2 (R"), 
and that 

dom (Ll' 2 ) = dom ((^:) 1/2 ) = Vl^ 2 (R"), (3.42) 

and hence provides precisely (3.17) in the particular case fl = R™, n G N. Quanti- 
tatively, one has 

\\ L a 2 f\\ L ^) ~ l|V/]|^(R»)n, (3.43) 

where the implicit constants depend only upon ai, a2, and n (and we recall that 
the symbol « indicates equivalent norms). 

Thus, Lemma 3.7 and (3.42) together with Theorem 3.6, in the particular case 
fi = R n , n G N, imply Theorem 3.9 below. However, we prefer to provide a second 
and alternative proof, reducing it directly to Corollary 2.6 as follows: 

Theorem 3.9. Assume Hypotheses 3.1 and 3.2 with fi = R" (so that W(fi) = 
W^' 2 (R™)), neN, and Zei L a denote the operator defined in (3.14). Then 

dom ((L a + q V) 1 ' 2 ) = dom (((i a + q V^)*) 1 / 2 ) 

v ; (3.44) 

= dom (Ly 2 ) = dom ((i:) 1/2 ) = iy 1 ' 2 (R"). 

Proof. Again, without loss of generality, we may assume in the following that Voo = 
and hence V 6 L P (R") for some p > re/2, p ^ 1. 

For the sake of notational convenience, we shall adapt the following convention 
for the duration of the proof of Theorem 3.9: The symbol e will denote a generic 
positive number (typically small, but not necessarily so), which may differ from one 
occurrence to the next, but which depends upon its value in the previous occurrence. 
We also frequently denote an identity or injection operator simply by the symbol I 
when the underlying spaces arc understood. 

Using this convention, we write V 6 L^ +£ (M. n ), and we factor V = V1V2, where 
Vj G L' i+£ (R n ), j = 1,2. Here, we have assumed that p is finite, as the case 
V G L°°(R n ) is simpler and is hence omitted in the following. We let B* and A 
denote, respectively, the operators defined by multiplication by V\ and V%. To prove 
the theorem, it suffices, in view of (3.42), to verify the hypotheses of Corollary 
2.6, with K = L 2 (R n ) = H, T = L a , and T = L a + q V. As we have noted 
above, Hypothesis 2.3 (i) is the principal theorem of [13]. To verify the remaining 
hypotheses, we begin by recording some simple observations. We set 

2* := 2n/(n - 2) (which equals 00 when n = 2), and 2* := 2n/(n + 2). (3.45) 

We then have by Holder's inequality, 

A : L 2 {R n ) -> L 2 * +£ (R n ), A : L 2 '- e {R n ) -> L 2 (R"), (3.46) 
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and similarly for A* , B and B* . Here, when n = 2, we interpret "2* — e" to be an 
appropriate finite exponent (to be precise, it is dual to the exponent 1 + e). Let 

L^(R n ) := {/ e L 2 (W l ) | (-A)"/ 2 / G L 2 (E")}, a > 0, (3.47) 

denote the usual fractional order Sobolcv space. For an appropriate choice of a G 
(0,1), we have, by (3.42), the Hardy-Littlcwood-Sobolcv theorem of fractional 
integration (see, e.g., [87, Theorem V.l], and (3.46), that 

dom (L)J 2 ) = W 1 ' 2 (R n ) C L 2 a (R n ) C L 2 '- e (R n ) n L 2 (R") C dom(A). (3.48) 

We obtain the same containment with B in place of A, and with L* in place of 
L a . Thus, (2.22) holds. Next, we verify (2.23) and (2.24). We recall that the heat 
semigroup e~ tLa enjoys the "hypercontractive" estimate (cf. [9, Corollary 3.5]), 

lle^'/IL^^^^'^ll/llLPCH-), 2*<p<g<2*. (3.49) 

Here we used the convention that < abbreviates ^ C for an appropriate constant 
C > 0. In this particular context, the underlying constant depends on p, q, n, a\ 
and a2, see [9]. (We remark that if n 3, then one can improve this interval by 
an e at both endpoints, but we do not require that fact here.) We observe that for 
t > 0, we have by the functional calculus that 

(L a + rl)^ 2 = jT ±- 2 e -C^ e-"-e-~. (3.50) 

Combining (3.49) and (3.50), we find that 

|| (L„ + T/) _1//2 || B ( L 2»+ e ( R „) ii 2( R „ ) ) + || (L a + r/) _1/2 || i3(i 2( R „) ii 2*- e ( R „)) 




(3.51) 



In turn, combining the latter bound with (3.46), and its equivalent for B* , and 
taking r = A + E, we obtain (2.23) and (2.24), with, say, R = 1 and E a = 0. This 
verifies Hypothesis 2.3 (ii). 

It remains to consider Hypothesis 2.1. Properties (2.1) and (2.2) are well-known. 
Property (2.3) follows from (3.46), and its equivalent for B* , and the fact that, in 
particular, the resolvent (L a — zl)^ 1 satisfies a hypercontractive estimate from 
L 2 * +eI (M™) to L 2 ~ e (W [ ) (see, e.g., [9]). Moreover, this hypercontractive bound 
holds with norm < \z\~ e for z G p(L) with — Rc(z) large, so that, for such z, the 
operator norm of K(z) is small; thus, property (2.4) holds as well. □ 

After focusing on the special case f2 = R™, n G N (which avoids the considera- 
tion of boundary conditions), we now turn to more general domains £1 C R n and 
to Dirichlet, Neumann, and mixed boundary conditions on dCl. To do so requires 
further assumptions on f2 and hence we will now assume that 51 is a strongly Lip- 
schitz domain. We recall that a strongly Lipschitz domain is a connected open set 
n C W" - (bounded or unbounded), whose boundary consists of finitely many parts, 
each of which is covered by a coordinate patch on which the boundary coincides (in 
an appropriately rotated coordinate system), with the graph of a Lipschitz func- 
tion, such that at most one of these parts is possibly unbounded. This includes, 
for instance, special Lipschitz domains (cf. [87, Sect. VI.3.2]), bounded Lipschitz 
domains, and exteriors of the latter (see also [19, Appendix B]). 
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Hypothesis 3.10. In addition to Hypothesis 3.1, suppose Q C R™ is a strongly 
Lipschitz domain. 

Remark 3.11. More generally, by virtue of the results in [26], one can consider 
domains which arise as bi-Lipschitz images of some smooth domain CI' C R™, where 
fi' may be R™, a half-space in R™, or a smooth bounded sub-domain of R n (or the 
complement, in either R" or a half space, of the closure of such). We refer the 
reader to [26] for details. In particular, it is possibly to replace the strong Lipschitz 
assumption in Hypothesis 3.10 by the assumption in this remark consistently in all 
subsequent results where Hypothesis 3.10 is employed. 

We note that the result of [13], the solution of the Kato square root problem 
for divergence form operators in R", has been extended to the setting of a strongly 
Lipschitz domain f2, where boundary conditions now play a non-trivial role. In the 
case that L a has Dirichlet or Neumann boundary condition, this was done in [25] , 
using a transference procedure from the global result of [13]. The case of mixed 
boundary condition was treated in [26] . More precisely, the result of [25] is that for 
f2 satisfying Hypothesis 3.10, Hypothesis 2.3 (i) holds with Tq = L 0) n,£> (Dirichlet 
condition), or Tq = L a ^.N (Neumann condition), H = L 2 (fl), and 

dom (Ll/lj,) = dom ((L^) 1 / 2 ) = << 2 (tt), 

(3.52) 

dom(L^ jiV ) = dom^L;^) 1 / 2 ) = W 1 ' 2 ^), 

where as usual Wq' 2 (Q.) is the closure of Cg°(Q) in W 1,2 (fi). The result of [26] is 
that 

dom {Ll% u ) = dom ((K^) 1 ' 2 ) = << 2 (tt), (3.53) 

where L a ,fi,n is defined in (3.11) and n C <90 is relatively open. 

We will now use the strategy behind Theorem 3.6 to extend Theorem 3.9 to 
the setting of a strongly Lipschitz domain and hence prove stability of square root 
domains in connection with Dirichlet, Neumann, and mixed boundary conditions 
on d£l. 

Theorem 3.12. Assume Hypotheses 3.2 and 3.10, let H C dQ be relatively open 
(n G {<9f2,0} permitted, cf. Remark 3.15 (ii)) and consider W '(O) = W^' 2 (f2). Then 
the following items (i) and (ii) hold: 

(i) V is infinitesimally form bounded with respect to -Ajj^, and there exist con- 
stants Mn tPt v,Q > and eo > (depending only on p, n and V and O) such that 
for all < e < Eo, 

\\\V\ 1/2 f\\l Hn) ^4(- A ^n) 1/2 f\\l m 

f 6 W^ 2 (Q). 

(ii) V is infinitesimally form bounded with respect to L a xiji and 

/ G Wn ,2 (n), < e < a^sa- (3.55) 
The form sum L a ji,n +q V is an m-sectorial operator which satisfies 
dom ((L„ A n +q V) 1 ' 2 ) = dom (((£„, n ,n +q V)*) 1 ' 2 ) 
= dom = dom ((K^/ 2 ) = << 2 (tt). 
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Proof. Without loss of generality we put the L°°-part Voo of V equal to zero. Since 
n is assumed to be a strongly Lipschitz domain, the Stein extension theorem (cf., 
e.g., [1, Theorem 5.24] or [87, Theorem 5 in §VI.3.1]) guarantees the existence of a 
total extension operator (cf., e.g., [1, Definition 5.17]) £ satisfying 

£ : W m > 2 {£l) W m > 2 (R n ) 7 



\\£ f\\ 2 wm , Hm ^ C m>n \\f\\* Q) , f G W m ' 2 (fl), m e {0, 1}, 



(3.57) 



(3.59) 



2(a)! 

for some constants C m .n > 0, m = 1, 2, and 

{£f)(x) = f(x) for a.e. ieS],/e L 2 {Vt). (3.58) 

Let V ex t and w ex t denote the extensions of V and v, respectively, to all of R™ 
defined by setting V cx t and v cx t identical to zero on R n \f2. Evidently, V ex t € L p (R n ), 
and consequently V cx t satisfies the hypotheses of Lemma 3.7. As a result, there 
exists a constant M n . P .v > for which 

|Kxt/||' 2(Rn) ^ £||(-A) 1/2 /||i2( Rn) + M n . p ye~ n/( - 2p ~ n ' , \\f\\ 2 L 2( Rn y 

f e W 1 - 2 ^ 1 ), e > 0. 
Consequently, using (3.57) and (3.58), one estimates 
ll«/lli»(n) = Kxt£/||| 2(K ») 

*S e ||(- A ) 1/2 ^/||i2 (R „) + Mn,p,ve~ n/{2p ~ n} \\£f\\ 2 L 2 (Rn} 
= eimiWlw + M n , P ,ve- n /^-^\\£f\\l 2m 
< eC 1; n||V/||| 2(n) „ + (eC^n + Co,oAf„, p ,ye-"/(^-"))||/||l 2(n) 
^ eCi,n||V/||| 2(fi) „ + (C 1J2 + C ,oA/„^y)6-»/( 2 f-»)||/||| 2(o) , (3.60) 

/ e w^ 2 (n), < e < 1. 

The form bound in (3.54) now follows by choosing e = eC^n throughout (3.60) and 
noting that 

l|V/||! 2(n) n = IK-A^n) 172 /!! 2 - w / e < 2 (^), (3.61) 

by the 2nd representation theorem (cf., e.g., [68, VI. 2. 23]), see (3.3) with a(-) = I n . 
Finally, we note that the results of [25] and [26] guarantee that (cf. (3.53)) 

dom (Z^ >n ) = dom ((L:,^ n ) 1/2 ) = < 2 (^)- (3-62) 
In light of (3.54) and (3.62), item (ii) now follows from Theorem 3.6. □ 

Remark 3.13. Without providing details, we note that in order to prove Theorem 
3.12 one can alternatively follow the proof of Theorem 3.9 essentially verbatim. 
The key observations are: first, that the containment (3.48) continues to hold 
since elements of W 1,2 (f2) may be extended to all of 1™ as discussed in (3.57), 
(3.58); second, that the hypercontractive estimate (3.49) continues to hold, with 
L a ,Q,D, L a ,(i,N ° r L a siji in place of the global operator L a , as does the analogous 
hypercontractive estimate for the resolvent. 

The following result improves on the infinitesimally bounded aspects in Theorem 
3.12 and proves the analog of input (3.18) in Theorem 3.6 assuming V satisfies the 
critical L p -indexp = n/2, n 3, in Hypothesis 3.3; at the same time it extends the 
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result of [33, Theorem 1.8.3] to the case of mixed boundary conditions on strongly 
Lipschitz domains: 

Theorem 3.14. Assume Hypotheses 3.1 and 3.3. Then the following items (i)- 
(iii) hold, and in each case r\ and rj are non-negative functions defined on (0, oo) 
possibly depending on Q, n, and the choice of boundary condition: 

(i) V is infinitesimally form bounded with respect to — Ajj.d, 

\W V \ 1/2 f\\ 2 mn) ^4(- A ^D) 1/2 f\\ 2 L2m +v(e)\\f\\h { n), f £ W^ 2 (n), e >0. 

(3.63) 

As a result, V is infinitesimally form bounded with respect to L a ^ o, 

hf\\l H n) < £Re[q 0)< , 2(n) (/,/)] +^( £ )||/||| 2(n)! / 6 1< 2 (11), e > 0. (3.64) 

(ii) If n = 1 (and hence SI is an interval), then V is infinitesimally form bounded 
with respect to — An,jv ; 

Hl^| 1/2 /||^(o) ^ e IK- A ^) 1/2 /|lL 2 (n)+^)ll/lli 2 (o), f&W 1 - 2 (n),e>0. 

(3.65) 

As a result, V is infinitesimally form bounded with respect to L a .n,iV; 

\\\ V \ 1/2 f\\ 2 L H n) < ^[q aiiyM(n) (/,/)] +»?(e)||/||£ a(n)) / G ^ 2 (0), £ > 0. 

(3.66) 

(in) Assume, in addition, Hypothesis 3.10, n 3, and let II C cT2 be relatively 
open (II G {<9S!,0} permitted, cf. Remark 3.15 (ii)). TTien V is infinitesimally form 
bounded with respect to — A^n, 

lll^| 1/2 /||Li (n) <e||(-An 1 n) 1/2 /||Li (n) +'?(e)||/|li a(n) , /G V^ 2 (fi), e>0. 

(3.67) 

^4s a result, V is infinitesimally form bounded with respect to L a jui, 

II l^| 1/2 /|lL ( o) ^ eRc Kw^ m (/./)] + *?( £ )H/llzW /e w£ 2 (fi), ^>0. 

(3.68) 

Proof. As before, we put the L°°-part VJx> of V equal to zero, 
(i). The infinitesimal form bound statement (3.63) may be found, for example, 
in [33, Theorem 1.8.4] in the case n 3 under the stronger assumption that Jl 
is connected. However, a careful look at the proof of [33, Theorem 1.8.4] and 
the results needed there, [33, Theorems 1.7.1, 1.7.6, and 1.8.3], reveals that the 
connectedness assumption on Q. is not necessary. Indeed, the principal elements of 
the proof of [33, Theorem 1.8.4] are the fact that W ' (0) may be embedded into 
L r (tt), 2* = 2n/(n-2) (cf., e.g., [34, Theorem V.3.6]), and that C§°(fi) is a form 
core for — Aq,d; these facts do not require connectedness of O. For the case n = 1, 
the form bound in (3.63) may be found, for example, in [68, p. 345-346]. The form 
bound (3.64) follows from the uniform ellipticity condition on a. In fact, from (3.7) 
and (3.63), one infers 

MWUn) < o^dM^w^n) (/>/)]+»?(*) Il/ll£»<n), / e W^ 2 (Q), e > 0. (3.69) 
Choosing e = die, £ > 0, in (3.69) yields (3.64). 
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(ii) . For the form bound given in (3.65), we refer once again to [68, p. 345-346]. 
Subsequently, (3.65) and (3.7) together imply 

(/, f)}+v(e)\\ff LHn) , f G W^ 2 (n), e > 0, (3.70) 

and (3.66) follows, choosing e = a\s. 

(iii) . We closely follow the proof for the corresponding Dirichlet case given in [33] 
and mentioned above in the proof of item (i). By a special case of the Sobolcv 
embedding theorem (cf., e.g., [1, Theorem 4.12, Part I, Case C]), 

W 1,2 (fi) ^ L 2 ' (ft), 2* = 2n/(n - 2), (3.71) 

where abbreviates continuous (and dense) embedding, and hence there exists a 
constant c > such that 

ll/lli»-(n) < c(||V/||i a(a) „ + \\f\\h m ), f e W^(Sl). (3.72) 
Using Holder's inequality, (3.72) implies 

(/, \W\}) LHn) < c||lU|| L „/ 2(n) (||V/||| 2(0) „ + ||/||| 2(0) ), 

/ e W 1>2 (fi), W e L n ' 2 {tt). 

Next, let e > be given. Since V & L™ /2 (ft), there exist functions V n/2 . e G L" /2 (0) 
and Uoo, £ 6 L°°(ft) with 

l|K»/2,e||i«/3(n) < e / c > ^0*0 = V n/2 , E (x) + Voo^{x) for a.e. (3.74) 
Applying (3.73) with W = V n / 2 .E, one estimates 

HIU'(n) = (/, \V\f)mn) < (/, [|K/2, e | + ||Voo ie || £ cc (n) ]/) £acn) 

<£||V/|| 2 iW +r/( £ )||/|| i2( n), (3-75) 

with 

r?(e) :=e+||Voo ie ||L-(n). (3-76) 

Noting the fact that 

l|V/||! w = ||(-An,n) 1/2 /|| i2(n) , / G Wft a (fi), (3.77) 

by the 2nd representation theorem (cf., e.g., [68, Theorem VI. 2. 23]), and using that 
Wji (Q) is a closed subspace of W 1,2 (f2), (3.75) then also yields 

Hlli» ( n) < e||V/||ia (n) - +»7(e)||/||ia ( n). / G w£ 2 (Q). (3-78) 

Since e > was arbitrary, (3.67) follows. 

To prove (3.68), one notes that the uniform ellipticity condition on a implies 

II V/|| iW WRe[q o ^, 2(n) (/,/)], /eW£' 2 (n). (3.79) 

Taking (3.79) together with (3.67) and (3.77), one infers that 

KH^ ( n) < or 1 cRfi[q. lW ri- a (n)(/'/)] + '?( c )ll/lli a (n). / e <' 2 (^). e > 0- (3.80) 
The form bound in (3.68) follows by taking e = oie, e > 0, in (3.80). □ 
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Remark 3.15. (i) The main idea for the proof of Theorem 3.14, is taken from the 

proof of the corresponding Dirichlet case in [33, Theorem 1.8.3]. 

(ii) Due to our strongly Lipschitz hypothesis on Q in Theorems 3.12 and 3.14 (Hi), 

W^(n) = W 1 ' 2 ^), (3.81) 

as noted in Remark 3.5 (Hi). Consequently, taking II = <9f2 in Theorems 3.14 and 
3.12 shows that (3.65) and (3.54) (resp., (3.55) and (3.56)) hold for the Neumann- 
type Laplacian — An,jv (resp., L a ^.N)- Alternatively, taking II — yields 

W^ 2 (Sl) = <' 2 (tt), (3.82) 

and shows that (3.65) and (3.54) (resp., (3.55) and (3.56)) hold for the Dirichlet- 
type Laplacian -An,n (resp., L a ,n,£>)- 

(Hi) Estimates of the type (3.63)- (3.68) can also be used to include first-order 
terms and hence consider the operator — div(aV • ) + (Bi ■ V ■ ) + div (B 2 ■ ) + V 
on n, assuming Bj G L n+s (Sl) n + L°°(f2) n , j = 1,2, for some e > 0. Rather 
than repeating the analysis for this more general situation, we refer to our detailed 
treatment of the critical L p -index for the case f2 = K™ in Section 4 (see also item 
(iv) below). 

(iv) In contrast to Theorem 3.12, where the analogs of 77(e) and 77(e) are explicitly 
provided in terms of Ce~ n / Iy2p ~ n \ the extension to the critical L p -index p = n/2 
in Theorem 3.14 no longer permits one to determine explicit estimates for 77(e) and 
77(e). Consequently, we can no longer use (3.20) in Theorem 3.6 (ii) to conclude the 
stability of square root domains as in (3.56). A systematic extension of this circle 
of ideas to the square root problem for the more general situation 

- div(aV • ) + (§! ■ V • ) + div (B 2 ■ ) + V on R" (3.83) 

at the critical L p -index will be undertaken in the following Section 4. 

(v) Introducing the set 

2>n(n) = I / G C °°(M"); supp (/) n (dQ\U) = 0} c W 1 ' 2 ^), (3.84) 

one recalls that || ■ ||\y I . 2 (n) an d j| • || q wn are equivalent norms on W(f2) (cf. 

Remark 3.4(H)), and hence on W^ 2 (fl). Thus, one infers that 

2?n(^) is a form core for — An,n and £a,fi,rij and hence also for £ a ,Q.n +q V, 

(3.85) 

employing the infinitesimal form boundedness in (3.68) (see also [34, Theorem 
IV.5.1], [68, Theorems VI.1.21, VI.1.33]). 

4. The square root problem for - div(aV • ) + (Bi • V ■ ) + div (B 2 ■ ) + V 

AT THE CRITICAL INDEX 

In this section we extend Theorem 3.9 to include first-order differential operator 
terms of the type (B\ • V • ) + div (B2 ■ ) in addition to the zero-order potential term 
V and we impose the critical (i.e., optimal) L p -conditions on Bj, j — 1, 2, and V. 

For this purpose we introduce the following hypotheses: 

Hypothesis 4.1. In addition to Hypothesis 3.1 with Q = W n , n G N, n ^ 2, 
suppose that 

B 1 ,B 2 eL n (R n ) n + L° c (R n ) n , V eL n / 2 (M. n ) + L° c (M. n ). (4.1) 
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Setting again L a := — div(aV • ), the operator formally given in (3.41) and pre- 
cisely defined by (3.14), we now consider the operator 

M a := L a +, (§! ■ V • ) + q div (B 2 • ) + q V, (4.2) 

in i 2 (R"). 

Decomposing Bj, j = 1, 2, and V as in (4.82) proves that under hypothesis (4.1), 
(Bi ■ V ■ ), div (i?2 • ), and V are all infinitesimally form bounded with respect to 
L a , rendering the form sum M a in (4.2) well-defined. 

For simplicity of notation, we will henceforth assume without loss of generality 
that the L°°-parts of Bj, j = 1, 2, and V in (4.1) equal zero. (It will be clear from 
the proof of Theorem 4.5 that these L°°-parts can be handled in a straightforward 
manner.) 

We start with the following preliminary result, in which we suppose that B\, B 2 , 
and V are small in the corresponding L p -norms, as in (4.3) below: 

Theorem 4.2. Assume Hypothesis 4.1 with £1 = W l , n ^ 3, and let L a and M a be 

defined as in (4.2) on W n , with Bj e L n (R n ), j = 1,2, V € L n / 2 (M. n ), supposing 
in addition the smallness condition 

ll^l|L»(R»)- + II^IL»(*»)- + WVh^^ < e . (4.3) 
If < £o = £q(u, ai, a 2 ) is sufficiently small, then 

dom(A/ a 1 / 2 ) = W 1,2 (R n ). (4.4) 

Moreover, 

\\Ma /2 f\\L^^)^C\\Vf\\ L 2 {Rn)n , f e W l ' 2 (R n ), (4.5) 
where C = C(n, ai ,a 2 , \\Bi\\ Ln{v , n)n , H^H^^^, II^IL"^^)) . 

Finally, these results continue to hold when n — 2, provided that Bj, j = 1,2, 
are divergence free, and that V = 0. 

Later in this section we will remove the smallness condition (4.3) (cf. Theorem 
4.5). 

Remark 4.3. (i) One observes that if B is divergence free, then div(Bf) = B ■ V/. 
Thus, in the case n = 2 it is equivalent to consider operators of the form M a = 
L a + B ■ V, with ||jB|| in(E n)« ^ £ o- 

(ii) Set R := M a - L a = (Bi ■ V • ) + (div B 2 -) + V. Then by Holder's inequality, 
(4.3), and Sobolev embedding, we have for n 3, with q fl (-, •) the scsquilinear form 
associated with R, and the decomposition V = uv, 

|q«(/,/)l = |(Bi7,V/) i2(B „ )n - (Vf,B 2 f) r2(Rnr + (u*f,vf) L2{Mn) \ 

<^o||V/||i 2(Rn) „. 

The analogous bound continues to hold for n = 2, but is slightly more delicate. 
In that case, by Remark 4.3 (i), we have that Rf = B ■ V/, which belongs to the 
Hardy space i? 1 (R 2 ), since B is divergence free. Indeed, by the result of [28], and 
(4.3), we have 

\\B-Vf\\ H1(m2) <Ce \\Vf\\ L z m 2. (4.7) 
Moreover, the endpoint Sobolev embedding result in M 2 yields that 

\\f\\BMO(W) < C||V/|| L 2 (R2)2 . (4.8) 
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Thus, by Fefferman's duality theorem [38], one has 

= \(B*f,Vf) L2{R2)2 \ <eo||V/||| a(B3)a . (4.9) 

Proof of Theorem 4.2. The proof will utilize the technology of the solution of the 
square root problem for L a developed in [13], [56], and [58]. We begin by observing 
that for £o sufficiently small, the operator M a and its Hcrmitian adjoint M* are 
m-accretive, by Remark 4.3 (ii) and the theory of sectorial forms and m-sectorial 
operators (cf. [68, Theorem V. 3. 2, Theorem VI. 1.33, Theorem VI. 2.1]). Moreover, 
setting 2* := 2n/(n + 2), 2* := 2n/(n — 2) as before, and more generally, setting 
:= pn/(n + p), p* := pn/(n — p), and p** := pn/(n — 2p), we arrive at the 
following result. (For simplicity, we abbreviate Il 2 (R") by I in this section.) 

Lemma 4.4. For < £o = £o( n , o,\, 02) sufficiently small in (4.3), the resolvent 
operator (I + t 2 M a )~ l enjoys the following properties. First, if n ^ 3, then 

|| (J + t 2 M a ) _1 || B( w)| w)) < C p Vp-i/2) , 2,- £l <^2, (4.10) 

for some < ei = £i(n, oi, a^)- When n = 2, £/ie previous bound continues to hold 
with 2* ^ p ^ 2, where the endpoint case p = 2* = 1 is interpreted as 

||(/ + t 2 Af a ) 1 || ;B ( ir i( K 2) ii 2( K 2)) ^ Ct ■ (4-11) 

In addition, for n ^ 2, 

\\tW + t 2 M a )-% (L2(mtL2(Mn)n) <C, (4.12) 

and /or n ^ 3, 

||t(J + t 2 ^)- 1 div || B(£ , (Rf . ) „ iia{aB)) < C q r"(V«-V2) , g = 2 - £2 , (4.13) 

for some < 62 = £2(h, ai, a 2 ). Moreover, for every pair of cubes Q\^Qi C W 1 , 
with side lengths £(Qi),l(Q2) satisfying 

0<t^ M[Qi), j = 1,2, and max(£{Q 1 ),£(Q 2 ) «S 100 dist(Qi, Q 2 ), (4.14) 

for n 3 we have the "Gaffney-type" estimate 

W + t Ma ) Hb(£ 2 (<2i),£ 2 (Q 2 )) + + * M a ) 1 || B(L 2 (Ql) i 2 ( Q 2) „ ) 

+ ||t(/ + i 2 M a )- 1 div|| B(L2(Qi) „ L2(Q2)) 

^ C , e -cdist(Q 1 ,Q 2 )/t ) k e R ( 4 15 ) 

For n = 2, 

i2iif 1 IUr7/r 1 j.2 n ,r \ — II 



|(7 + t M a ) || B ( i 2(Q 1 ) )i 2(Q 2 )) + pV(/ + t M a ) || e ( L 2(Q 1 ) i 2(Q 2 )2) 

+ \\t(I + t 2 Ma)- 1 div 



B(L 2 (Q 1 ) 2 ,L 2 (Q 2 )) 

< C ( T- 7F> — 7T\ ) » fc e N - ( 4 - 16 ) 
Vdist(Qi.Q 2 )/ 

T7ie various constants c, C, C p tmd C q in estimates (4. 10) -(4. 16) depend upon n, 
Oi, afid a2- 

We remark that there is nothing special about the 4th order decay in the two- 
dimensional Gaffney estimate (4.16): indeed, in the proof of this part of Lemma 
4.4, we could obtain polynomial decay of any specified order V, by choice of e 
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small enough depending on N, but as it turns out, 4th order decay will suffice for 
our purposes here in two dimensions. 

Let us take the Lemma 4.4 for granted momentarily. We shall defer its proof to 
the end of this section. Before proceeding further, let us note for future reference 
the following immediate consequences of (4.15). First, by (4.15) and Sobolev's 
inequality, we have for n 3 that 

11(7 + t 2 M < +-l p -cdist(Qi,Q 3 )/t (4 171 

where c and the implicit constant depend on n, a±, aa, ||-Bi||l"(K")™, ||-S2||l™(r™)™i 
and ||V"|| i „/2 (M „ ) . Next, given a cube Q C R™, set S k (Q) := 2 k+1 Q\2 k Q, and let 
< t ^ 4£(Q). We then have the following variant of the "Gaffncy" estimate, for 
n 3: 

\\(l + t M a ) || B (£2(5 fc (Q)) )i 2(Q)) + ||fV(J + I M ) || B (£2( Sfc (Q)) ii 2(Q)n) 

+ ||i(J + t 2 M a ) 1 div || e(i 2 (Sfc( Q )) „ ii 2 (Q)) 
^ Ce-^^W)/*, fc e N, (4.18) 

as one may deduce from (4.15) by covering Su{Q) by a bounded number of cubes 
Q' of side length £(Q') w 2 fc €(Q), with dist(Q, Q') w 2 k £{Q). Similarly, from (4.17) 
we obtain for n ^ 3 that 

\\(l + t 2 M 1 _1 II < +-i - C 2^(o))/t c 419 i 

When n = 2, we have 

+2 JIT \ — 1 



(I + t Af ) || Z 5 (L 2 (Sfc( Q ))i 2 ( Q )) + ||lV(7 + I A/ a ) || e ( L 2 (Sfc( Q ))i 2 ( Q ) 2 ) 



+ \\t(I + t 2 M a )~ 1 div\ 



t ^ 



* c {*m) ' fceN - (4 - 20) 

1 /2 

Returning to the proof of Theorem 4.2, we resolve M a as follows: 

/>OG 

M^ 2 = c / di (J + t 2 M a )- 3 t^Mlt' 1 (4.21) 
./o 

We set F a := L a + (£?i • V • ), so that, for n ^ 3, we have M = H a + div(B 2 -) + V, 
and in two dimensions we simply have M a = H a . Thus, for n 3, we have, 
1*00 

Ml' 2 = c dt(I + t 2 M a y 3 t z M a H a t~ x + M- 1 ' 2 [div(B 2 • ) + := Ti + T 2 , 
Jo 

(4.22) 

1/2 

and when 77. = 2, we have M a = Ti; the term T 2 is vacuous. Our goal is to show, 
in all dimensions n J? 2, that 

llTiFH^^) ^C||VF|| L 2 (Rn) „, Fe^fl"), (4.23) 

and for 77 ^ 3, that 

llTaFllia^^CeollMy^l^llVFll^^n, F e W 1,2 (R n ). (4.24) 

Here we abbreviated, 

ll M a 1/2 || op := sup ||mV^|| (4.25) 
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Under the a priori assumption that \\M a \\ is finite, the bound (4.5), and there- 
fore also the conclusion of Theorem 4.2, then follow for Eq small enough. One 
may remove the assumption of a priori finiteness by temporarily restricting F 
to the class Cq°(B(0, R)), and temporarily replacing Bj, j = 1,2, and V by 
Bj,N '■— BjXng i<jv}' 3 ~ an< ^ '■— min(V,iV), under which circumstances 
we have a priori finiteness by [23, Proposition 11, p. 19]. We then establish (4.5), 

1/2 

with bound independent of R and N, for the associated M N , at which point one 
may pass to a limit. We omit the details. 

We now turn to the proofs of (4.23) and (4.24). 

The bound (4.24), which arises only when n 3, is the easier of the two, and we 
treat it first. Recalling that 2, := 2n/(n + 2), 2* := 2n/(n — 2), we record several 
simple observations: 

||F|| i2 . (a „) ^C\\VF\\ L 2 (Rn) n, (4.26) 
\\VF\\ L2 , iRn) ^e \\F\\ L2 * {mn) , (4.27) 

\\Bi ■ ,f\\ L 2, {Rn) < eo||/1| L 2 (R „ ) „, (4-28) 

\\^F\\ L2{Rn)n <e \\F\\ L2 * (Rn) . (4.29) 

The first of these is just Sobolev embedding, and the other three follow directly 
from Holder's inequality, given (4.3). Moreover, we claim that 

ll M a" 1/2 /|L 2(R „) < C||^a 1/2 || op ll/ll^(R«)> (4-30) 

||iV^- 1 / 2 div(/)|| £3(R „ ) CllM^H^II/lU^n. (4.31) 

Taking this claim for granted momentarily, one may then obtain (4.24) by combining 
(4.26), (4.27), (4.30), and (4.31). We omit the routine details. Let us now establish 
the claimed bounds (4.30) and (4.31). The former is equivalent to the dual estimate 

\\(M:r 1/2 f\\ L2 * (Rn) < CllM^IIJI/ll^), (4.32) 

where M* is the adjoint of M a , while (4.31) is dual to 

\\y(M:)- x / 2 f\\ L ^ n)n < CllM^II^II/ll^n). (4.33) 

Moreover, (4.33) implies (4.32), by Sobolev embedding, so it suffices to verify (4.33). 
By cllipticity, 

||V(K)- 1/2 /|r i2(R „ )n < (a*V(M:)-^.f,V(M:)-^f) L2{Rn)n 
= q La {(K)- 1/2 f,(M* a )- 1/2 f) 

= ((M* a y/ 2 f,(M* a )-v 2 f) L2Q!in) - (B;(AP a )-^f,v(M:r^f) L2{Rn)n 
(v.*(Kr 1/2 fMKr 1/a f)vw + (v(M:)- 1/2 /,s 2 (A/:)- i / 2/ ) L2(Rn)n 

=: I + 11 + III + IV. (4.34) 

For £o small enough, we may use (4.27), (4.28), (4.29), and Sobolev embedding to 
hide // + III on the left-hand side of the inequality. For the main term, we have 

|/| = |(/,MV 2 (K)-vv) L2(Kn) | < ll/ll^jll^iy v(m;)-vv|| £W ., 

(4.35) 

whence (4.33) follows. We have therefore established (4.24). 
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It remains to prove (4.23). It is enough to show that for g, F <G Cq°(R"), we 
have the estimate 



1(5,7^)^)1 < C\\g\\ L , {Rn) \\VF\\ L , (R n )n . 



(4.36) 



By the definition of T\, and by Cauchy-Schwarz, it is enough to establish the square 
function estimates 



and, for n ^ 3, 



t a M* I + t 2 M* 



t(l + t 2 M a ) 1 H a F 



9 



2 d n xdt 
t 



< C\\g 



li 2 ( 



2 d n xdt 
t 



< C||VF||i a(RB) , 



(4.37) 



(4.38) 



In two dimensions, we simply change H a to M a in (4.38). Since M a is a one-to-one 
maximal accretive operator (given the small norm condition (4.3)), it therefore has 
a bounded holomorphic functional calculus (sec [2, Theorem G, p. 102]), whence 
estimate (4.37) follows directly from the Mcintosh- Yagi Theorem (see, e.g., [2, 
Theorem F, pp. 98-99] or [76].) We therefore turn to the proof of (4.38), which is 
the heart of the matter. Let us suppose for now that n ^ 3; we shall discuss the 
modifications required to treat the two-dimensional case at the end of the proof. 
The first step is the standard reduction to a "Carlcson measure" estimate, following 
the analogous step, with some modifications owing to the presence of lower order 
terms, in the treatment of the square root problem for L a in [13], [14], [23], [56], 
[58]. 

Before proceeding, we introduce two standard averaging operators that will be 
used in the sequel. Let P t be a "nice" approximate identity, that is, 



Ptf 



*f: 



(4.39) 



with 4>t(x) := t~ n <fi(x/t), for some non-negative, radial <f> <E Cq°(R™), supported in 
the unit ball in R™, with integral equal to one. Let E t denote the dyadic averaging 
operator, that is, if Q(x,t) is the minimal dyadic cube (with respect to the grid 
induced by Q) containing x, with side length at least t, then 



(Etg)(x) := 



1 



d n yg(y)- 



\Q(x,t)\ J Q ( X:t ) 

Returning to the proof, we let 9 t denote the operator defined by 

-l 



t 2 M a ) 



(div(o) - B 1 ) ■ f, 



t 2 M, 



Otf ■■= -t {I 

and observe that therefore 

t(l + 

Using a strategy that originates in 
t(l- 

= (6 t - 9 t I n {x) ■ E t P t )WF{x) + 9 t I n {x) ■ E t P t {VF)(x) 
=: R t {VF){x) + 9 t I n (x) ■ E t P t {VF){x) =: I + II, 



,) H a = t V . 

, we may then write 



(4.40) 

(4.41) 
(4.42) 



t 2 M a ) l H a F{x) = 6 t VF{x) 



(4.43) 
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where I n is the nxn identity matrix, and where Et and Pt have been defined above. 
The key term here is II. By Carleson's embedding lemma, the desired bound for 
II follows immediately from the Carleson measure estimate 

sup ±£ \(6 t I n )(x)f^^C, (AAA) 

where the supremum runs over all cubes Q C R n , and where Rq := Q x (0, £(Q)) 
is the standard "Carleson box" above Q. To establish (4.44) is our main task, and 
we postpone the proof momentarily. 

Turning to the term J, we observe that by a standard "Schur's lemma" argument, 
and standard Littlewood-Paley theory for the classical heat semigroup, to obtain 
the desired bound for R t , that is, to show that 

t Ji? t Vif ^<C||VF||| 2(Rn) „, (4.45) 
*+ 

it suffices to establish the quasi-orthogonality estimate 

||i? t V S 2 A e 2s2A F|| i2(R „ )5 :Cmin^)",Q^ \\s e^ A V ■ VF|| L2(R „ ); (4.46) 

for some fixed /3 > 0, where again, A = J3" =1 d%. is the usual Laplacian on W 1 . 
To this end, let us observe that by (4.18), (4.19), and (4.28), we have that 9 t also 
satisfies estimate (4.18) (with slightly different constants), that is, for < t ^ A£(Q), 
we have 

\\0t\\ B{ LHs k (Q))",L* m <Ce~ c2H ^l\ ken. (4.47) 
Similarly, by (4.10) and (4.28), we have 

sup ||6» t || S ( L 2 (Rre) „ :L2(R „ )) sC C . (4.48) 
*>o 

In particular then, by [4, Lemma 3.11] we infer 

sup||(6» t J n )i; t || B(L 2 (R „ ) „ ii 2 (ffi „ )) <C. (4.49) 
t>o 

Moreover, (8 t I n )E t Pt trivially satisfies the same bound as t in (4.47), since the 
kernels of E t and Pt have compact support in a ball of radius t; thus R t also satisfies 
(4.47) and (4.48). Note further that by construction, Rtl n = 0. Consequently, by 
[4, Lemma 3.5(i)], 

||i? t V S 2 Ae 2i ' 2A F|| < Ct||V 2 ( S 2 Ae 2s2A F)|| L 

II Hi ( r ) II V /Hi/ (r ) 

= C(t/ S )||( S 2 V 2 e s - A ) se s - A A J F|| i2(R „ ) < C(t/ s )\\ S e s - A AF\\ L2(Rny 

which establishes (4.46) in the case t ^ s. 

Next, we verify (4.46) in the case s < t. Recalling that Rt is a difference, 

Rt = 6 t - (6tIn)E t Pt, (4.51) 

we now consider the contributions of the two terms separately. We have 

\\(9I n )E t P t V S 2 Ae 2s2A F\\ L2{Rn) = l\\(0I n )E t tVP t e s * A se a ' A AF\\ La ^„ y (4.52) 
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whence (4.46) follows for the present term, with (3 = 1, by (4.49) and the uniform 
I/ 2 -boundedncss of tVPt and of e s A . Consider now the contribution of 9 t . To 
simplify notation, we set Q s := s 2 Ae 2s A , and note that 

6 t VQ s F = -t (l + t 2 M a y l (dWia^-Bi^VQsF 

= t(l + t 2 M a y 1 M a Q s F -t^ + tHla)' 1 [div(B 2 -) + V] Q S F 
=: - IF , . (4.53) 

Then 

G s>t = -(/-(/ + t 2 M a )- 1 ) se 2s2A AF, (4.54) 

which, by the case p — 2 of (4.10), satisfies (4.46) with /3 = 1. To conclude our 
treatment of Rt, it remains to consider H s ^, which we split into 

H s . t = (I + t 2 M a ) 1 div (_B 2 Q S F) + (l + t 2 M a y X VQ s F =: H' Bjt + H'^ t . (4.55) 

By (4.10), with p = 2* — e, for some e > 0, and by Holder's inequality, we have 

\\Kth 2 m £ t 1 - n( - 1/p - 1/2) \\VQ t F\\z, 0l » ) < t^/v-^WQsFW^^ 

= t l-n(l/ P -l/2)UsiA sAe s>A F \ 



ILp**(R") 

< t l-«(l/p-l/2) s l-™(l/2-l/p")|| sAe;S 2 A^| 



L 2 (R") 



-1- 



sAe s2A F|L 2r v (4.56) 



where we recall that p** := pn/(n — 2p), and we have used the well-known bound 
||e s A ^|| ir ( R „-) ^ Cr s~ n ^^ 2 ~ 1 ' r > ||/i||l 2 (r») I0r r ^ 2; the last step is an elementary 
computation, using the relationship between p and p** . We observe now that 

n n (1 2 + n\ (\ 1 \ , , 

since p = 2* — e. Similarly, by (4.13), we have that for some q = 2 — e', 
KtlU'CR-) S t-" (1/9 " 1/2) \\B 2 Q S F\\ rq(B ^ < f- n (V«-Va) \\Q.F\\ Lt . 



_ .-n(l/?-l/2) II s 2 A S A„^ 2 A F || 

< f -n{l/q-l/2) l-n(l/2-l/q-*) II „A„s 2 A p || 



n(l/ 9 -l/2) , 

sAe F\\ T 2(mn\j (4-58) 



Since q <2, this concludes the proof of (4.46), and hence also that of (4.45). 

Modulo Lemma 4.4, it therefore remains only to establish the Carleson measure 
estimate (4.44). To this end, we invoke a key fact in the proof of the Kato conjecture 
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for L a - By [13], there exists, for each Q, a mapping Fq : W™ — > C n such that 
(i) f d n x\(VF Q )(x)\ 2 ^C\Q\ 



(ii) I d n x\(L a F Q )(x)\ 2 <C 



£(Q) 2 

i r e{Q) r „od n xdt ( 4 - 59 ) 



(m) sup^ / / |C'(M)| 2 
Q 



\Q\Jo J Q ' t 

sjCsup — / / IC^^V^qM 



q IQI Jo Jq t 

for every function £ : R™ +1 — > C™, where VFq := (1?^- (Fq) k)i^j,k^n is the Jacobian 
matrix, and where the product 

n 

(E t VF Q = CjE t DjF Q (4.60) 

is a vector. Moreover, although not stated explicitly in [13], the mapping Fq also 
satisfies a slightly stronger version of (4.59) (i), namely that for some q = 2 + e > 2, 
with e depending only on dimension and cllipticity, we have 

d n x\(VF Q )(x)\ q <tf|Q| (4.61) 

(see, e.g., [58] for a proof of this fact in the case that the kernel of the heat semigroup 
e~ tL enjoys a pointwise Gaussian upper bound, but the proof there may be readily 
adapted to the general case in which the pointwise kernel bound is replaced by one 
of "Gaffney-type" ) . 

Following the strategy in [23, Chapter 3], and given the existence of a family 
of mappings Fq with these properties, we see that by (4.59) (Hi), applied with 
C(x,t) = 8tl n , in order to prove (4.44), it is enough to show that 

*(Q) r d n rdt 

' (6 t I n ) (x) (E t VF Q ) (.t)| 2 ^— < C\Q\ , (4.62) 

uniformly in Q. To this end, we write 
(e t I n )(E t VF Q ) 

= (0 t I n ) ((E t - E t P t )VF Q ) + ((9 t I n ) EtP t - t ) VF Q + 6 t VF Q 

= : +R[ 2) +e t VF Q , (4.63) 

and notice that — R t = Rt, as defined in (4.43). Thus, the desired bound for this 
term follows immediately from (4.45) (with Fq in place of F) and (4.59) (i). We 
claim that a similar bound holds for the contribution of R^\ by (4.59) (i) and the 
well-known square function estimate 

d n Tf\f 

^ \(E t P t )h(x)\ 2 — r < C\\h\\^ n) . (4.64) 

+ 

Indeed, since E t is a projection operator, we have that 

E t -E t P t = E t (E t -P t ). (4.65) 
Consequently, we may invoke (4.49) to reduce matters to (4.64), with h = VFq. 
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The remaining term is 
6 t VF Q =t(l + t 2 AI a y 1 (L a F Q )+t(l + t 2 M a y 1 (B 1 -VF Q ) =: $ t + ¥ t . (4.66) 
By the case p — 2 of (4.10), and by (4.59) (ii), we have that 



rl n rrlf f * WJ f 

1 Jo Jts." 



i 

<<7|Ql w/o ^ = C|QI ' (4 ' 67) 

as desired. Finally, to treat the contribution of ^t, for = 2 + e as in (4.61), we 
choose p so that p* = pn/ (n — p) = <z, and note that since q > 2, we therefore have 
p > 2*. Then by (4.10), followed by Holder's inequality, 

rKQ) r fin Tf u rt(Q)/ r - \ 2/p 

/ \y t ( x )\*±J?E. < ( I^-VFqIVx) ti-anCVp-i/aJdt 

< / Hall nvf>ji 2 , ^-2«(i/p-i/2)j f 

— / H i II L^fR 71 )™ H Q II Lp (R™) " r 



/••«(Q) 

< \Q\ 2 'p" / t x-an(i/p-i/a) {ft) (4 68) 



where in the last step we have used (4.61) with q = p* . Since p > 2» = 2n/(n + 2), 
we infer 

1 - 2n(l/p - 1/2) > 1 -2n ( - M = 1 - (n + 2 - n) = -1. (4.69) 

\ 2rt 2 J 

Consequently, the last term in (4.68) equals 

C |Q|2/p*+2/n-2(l/p-l/2) = C |Q|. (4 . 70 ) 

This finishes the proof of Theorem 4.2, in the case n ^ 3, modulo the proof of 
Lemma 4.4. Before proving the latter, let us now treat the two-dimensional case. 

In two dimensions, it remains to verify the square function bound (4.38), but now 
with H a = M a = L a + B ■ (V • ), where div (i?) = 0. As in the higher-dimensional 
case, our first step will be to reduce matters to the Carleson measure estimate 
(4.44). The second step (again the main step), will be to prove that (4.44) holds. 
We define 9t as above (4.41), with B\ = B, but observe that now, 

t V = t(I + t 2 M a yHl a = t- 1 (/-(/ + i 2 M )- x ) =: t^Tt (4.71) 

(cf. (4.42)). We note that by Lemma 4.4, the operator T t satisfies 

sup||TJ e(L2(R2)) sCC, (4.72) 
t>o 

and also estimate (4.20), that is, for < t Ai(Q), we have 

||Tt|| B( L W Q)),L=(Q)) < C (afc^Qj) > fceN - ( 4 - 73 ) 

Let us now proceed with 

Step 1: reduction to the Carleson measure estimate (4.44). 
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Verification of Step 1. The goal here is to show that the square function bound 
(4.38), in two dimensions, follows from (4.44). To this end, we make the same 
splitting as in (4.43), recalling that in the present setting, H a = M a . Then main 
term II is handled exactly as before, and is bounded, without further ado, in terms 
of the Carleson measure estimate (4.44). The "remainder" term / is also handled 
as before: we need to verify the quasi-orthogonality bound (4.46). To carry out the 
latter task will require a bit of an argument. 
In the present setting, by (4.71), we have 

R t V = t V - (6 t I 2 )E t P t V = t^Tt - t-^Tt^EtPtV , (4.74) 

where I 2 = V<1> is the 2x2 identity matrix. We first establish (4.46) in the case 
s ^ t, and consider the two terms in (4.74) separately. The contribution of the 
term t -1 Y t = 6* t V is treated exactly like G s .t in (4.53); there is no term H s t in the 
present context. The contribution of the other term {9 t h)EtPt^ = t^ 1 (T t (f>)_B t P 4 V 
may be handled exactly like its higher-dimensional analogue in (4.52), except for 
one minor technical issue which we now describe. In estimating (4.52), we used the 
bound (4.49), which in turn relied upon (4.47) and (4.48), but in two dimensions it 
is not clear that 9 t satisfies the latter pair of bounds. Thus, we cannot simply invoke 
[4, Lemma 3.11] "off-the-shelf to obtain (4.49); nonetheless, the proof of [4, Lemma 
3.11] goes through, mutatis mutandi, to establish (4.49), since 9 t I n = i _1 T t <l> (cf. 
(4.71)), where T t satisfies (4.72), (4.73), and T t l = 0. The rest of the argument 
may now be carried out as before. We omit the details. 

Next, we consider the case t < s. As above, let Q s = s 2 Ae 2s A . We now claim 
that for g <G W 2 ' 2 (R 2 ), we have 

\\RtVg\\ LHM *) < Ct\\V 2 g\\ L2(R2y (4.75) 

Let us take this claim for granted momentarily, and apply (4.75) with g = Q S F. 
We then have 

II^V( S 2 Ae 2 ^F)|| L2(R2) < ^|| S 2 VV 2 * S e^AF|| i2(R2) 

* (4.76) 
< -\\se s2A AF\\ 

which is (4.46) in the case t < s, with ,3 = 1. 

We now turn to the proof of (4.75). For t fixed, choose j(t) so that 2^ f ' ^ 
t < 2 J (*) +1 , and let Dt denote the dyadic grid of cubes of side length 2^*>. For a 
measurable set E, let j E f := j E f denote the average of / over E. For (x, t) 

fixed, we set 

G x , t (z) := g(z) - g(x) -(z-x)- (V. 9 ) (x). (4.77) 
Observe that, by (4.74), 

R t Vg(y) = (T t - (T t $)E t P t V) (r 1 ^,*) (y) =: f t (t^G^) (y), (4.78) 

since T t l = 0, and since (i? t V)$ = 0. We note that T f satisfies (4.72) and (4.73) 
(for the term (T t ^)E t P t V , this relies on the fact that, as noted above, estimate 
(4.49) continues to hold in two dimensions). Thus, setting So(Q) := 2Q, and 
S k (Q) := 2 k+1 Q\2 k Q, we have 

\\RtVg\\ L 2 {m = (j2 [ \rt{t~ 1 G x ^(y)\ 2 d 2 y] / 
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1/2 



<£( £ 7 / IT^G^ls^K^dV 2 * 



Q-'Q 



1/2 



<£2"^ ^ / f d 2 yd 2 x 



\g(y)-g(x)-(y-x)-(Wg)(x)\ " 
2 2k t 



1/2 



(4.79) 



where we have used that £(Q) ~ t, for Q £ ©t, and that n = 2. By Plancherel's 
Theorem, the expression in curly brackets equals 

2 

d 2 hd 2 x 



h\<C2 k t 



\g{x + h)-g{x)~h-{Vg) {x)\ 



2 2k t 



, , A _ x _ {2m^-2 k th)\\ 2 l ,,, , Xl , , , 

'' LoU W ' ) 'f 1 '^'^ 

<t 2 ||A. g ||| 2(R2) . (4.80) 
Plugging this bound into (4.79), we obtain (4.75), which completes Step 1. □ 

Step 2: proof of the Carleson measure estimate (4.44). 

As in the higher-dimensional setting, we use the family of mappings Fq satisfying 
(4.59) and (4.61), so that once again matters are reduced to proving (4.62). We 
make the same splitting, in (4.63), and observe that the the "main term" O^Fq, as 
well as the two "remainder terms" iij and R[ 2 \ may be handled exactly as before. 
To handle the remainder terms, we use (4.49), which, as we have noted above, is no 
longer a direct consequence of Lemma 3.11 of [4] as stated, but nonetheless follows 
readily from its proof; we omit the details. Modulo the proof of Lemma 4.4, the 
proof of Theorem 4.2 is now complete. □ 

Next, we remove the smallness condition (4.3), and hence formulate the principal 
result of this section: 

Theorem 4.5. Assume Hypothesis 4.1 with n G N, n ^ 3 ; and let L a be defined 
as in (3.41), and M a be defined as in (4.2). Then 

dom(7\/ f y 2 ) = W 1 ' 2 ^). (4.81) 

Moreover, decomposing Bj, j = 1,2, and V into 

ll-Sj'jEOiOO || ^ooflgnjn < °°! j = 1)2, 

V = V e0tn /2 + V e0tOO , \\V £0tn/2 \\ L n/ 2{m) ^ £o/3, ||K ,oo|U-(R") < 00, (4.82) 
choosing < £o sufficiently small such that Theorem 4.2 applies to M a ^ EQ , where 
M aiS0 =L a + B hS0tn ■ V + div (B 2 , S0 ,n ■ ) + V eo>n/2 , (4.83) 
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and 

M a = M a>E0 + S 1j£0iOO ■ V + div (S 2i£0iOO ■ ) + V E0tOO . (4.84) 
Then for all e > 0, 

\\ M « /2 f\\Uu») < [Ci+eC 2 ]\\V.f\\ 2 LW + [C 2 /(4e)]||/|| 2 L2(R „ ) , / e WA 2 (]R"), 

(4.85) 

where 

Ci = Ci(n, ai, a 2 , ||-Bl,e ,n|| i n( K n)») ||- B 2,eo,n|| i „™„ :) „, ||K ,n/2||l,™/2(Rn)), 

(4.86) 

C 2 = C2(||-B liE0iOO || ioo(Ja „^ n , ||S 2ie0 ,oo|| iO c.( K „)„j \\V eo ,oo |U°°(R™))- 

Remark 4.6. We observe that, in Theorem 4.5, our quantitative bounds do not 
depend simply on ||Bj|| n and ||V|| n / 2 , but rather also upon the L°° norms of the 
"non-small" parts of Bj and V. Thus, we do not obtain a homogeneous estimate 
(i.e., one which involves only ||V/|| 2 ) on the right hand side of (4.85). One should 
contrast this situation with that of Theorem 4.2, in which we obtain the homoge- 
neous estimate (4.5), in the presence of (4.3). 

Proof. We note that the sequilinear forms uniquely associated with -Bi, Eo ,oo • V and 
div (J3 2 , eo ,oo • ) are of the type 

(^i,e ,oo/,V<?) L2(Rn) „, f,g€W^ 2 (W n ), (4.87) 

and 

-(V/ 1 ^2, eo ,ooff) £a(B „ ) » 1 /,ffGW 1,2 (R n ), (4-88) 
respectively. Thus, for b € L°°(K™), the elementary estimate 

\(bf,(-A)^ 2 f) L2{Rn) \ < 11611^^) [ £ ||(-A) 1 / 2 /|| 2 L2(R „ ) + (4e)- 1 |i/||i 2(K , 

/ g W rl - 2 (R"), e > 0, (4.89) 

shows that the forms (4.87) and (4.88) are infinitesimally bounded with respect to 
the form (V/, Vg)i2(i»)n , /, g € W 1 ' 2 (M. n ), and hence also with respect to the form 
uniquely associated with M a ^ Eo . 

At this point it remains to apply Theorem 2.14 repeatedly as follows: We identify 
T with -A defined on W 2 ' 2 (R n ). Then, in a first step, identifying A = V e0tOC , 
B* = l£2( R ™), T\ = M ajSo , and T = T\ + V E0lOO , it is obvious that Hypotheses 2.3 
and 2.13 are satisfied and hence Theorem 2.14 applies to T = M a , eo + V eo ,oo- In 
the second step we now identify A = V, B* = Bi tS0tOO , T\ = M a ^ + V e0tOO , and 
T = Ti+Bi tS0tOO • V, again Hypotheses 2.3 and 2.13 are satisfied and hence Theorem 
2.14 applies to T = M a<eo +V r £0iOO + -Bi. £0lO o -V. Indeed, the integrability conditions 
(2.23) and (2.24) are satisfied since there exists constants C\ > 0, C 2 > such that 

||(-A) 1 /2(_A + (A + i?)/ i2(R „ ;rf , l;c) )- 1 /2|| B(i2(Rn)) <Clj \ >RiE ^E Ql 

||pi,e ,oo|(— A + (A + i?)/ i 2 (R „ ;dre;E )) _1/2 || B ^ i2 ^ Rre ^ C 2 A~ 1/2 , A > R, E ^ E , 

(4.90) 

for < R, Eq sufficiently large. (Actually, to be precise, one has to apply this step 
by step to each term (Bi je0i0O ) .(V)j, 1 ^ j ^ n, but the necessary estimates are 

the same.) In the third and final step one identifies A = B2,e ,oo> B* = div(-), 
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Ti = Ma. So + K ,oo + (£i, Co ,oo • V • ) , and T = M a = Ti + div (S 2 , £0lOO • ) and once 
more uses the estimates in (4.90). □ 

Remark 4.7. (i) We note that the results of Theorem 4.5 continue to hold when 
n = 2, provided that Bj permit a decomposition as in (4.82) such that i?j. £o ,2, 
j = 1,2, are divergence free, and that V = 0. 

(ii) The critical condition V € L 1 (R 2 ) is not sufficient in the case n — 2 as the 
Sobolev embedding (3.71) fails in this case. Thus, for n = 2 one resorts to the 
condition V € L P (M. 2 ) for some p > 1 as in Hypothesis 3.2. 

(iii) While on the L p -scale the critical index for V in Theorems 4.2 and 4.5 equals 
p = n/2 for n ^ 3, it is possible to go beyond L n ^ 2 (fl). In this context we refer to 
Yagi [91] (see also [90]) who proved stability of square root domains in the particular 
case of three-dimensional Schrodingcr operators, that is, in the special case a = I3, 
and in the case Bj = 0, j = 1,2, for potentials V in the Rollnik class 72-3, that is, 
measurable functions V : R 3 — s- C satisfying ||T^||k 3 < 00, where 

\\V\\tz 3 := [ cPx<Px>\V(x)\\x-a/\- 2 \V&)\. (4.91) 

By Sobolev's inequality, L 3 / 2 (R 3 ) C 72-3, in fact, for some C > 0, 

\\V\\n 3 ^ C\\V\\ L3/HSS) (4.92) 

(cf. [85, Section 1.1]). In fact, Yagi [91] also considers fractional powers a € (0, 3/4), 
not necessarily a = 1/2, in this context. 

Finally, we turn to the proof of Lemma 4.4: 

Proof of Lemma 4.4. We begin by proving (4.10), (4.12) and (4.13). We shall ob- 
tain these estimates by establishing appropriate convergence of a Neumann series 
expansion of the resolvent. Suppose first that n ^ 3 and set 

R := (Bi • V • ) + div (B 2 -)+ V, (4.93) 

so that 

(I + tHla)- 1 (4.94) 
= (/ + t'La)- 1 / 2 [I + ((/ + t 2 L a )-^ 2 t 2 R (I + t 2 L a )-^ 2 )] -\l + t 2 L a )-^ 2 

OO 

= (J + t 2 L a )- 1 / 2 J2 [ ~ (I + t 2 L a )- l/2 t 2 R (I + t 2 L a )- 1 ' 2 ] \l + t 2 L a )- 1 / 2 . 

k=0 

We shall require the following estimates, which may be gleaned from the results in 
[9]: 

\\(I + t 2 L a )-^g\\ LqiMn) < C^i-^/p-V*) \\g\\ w , (4.95) 

2* -e sjp q ^p* ^ 2* +e, 
\\ty(I + t 2 L a )-^ 2 g\\ LPiRn)n ^C p \\g\\ LPiRn)l 2* -e^p^2 + e, (4.96) 

\\t(I + t 2 L a y 1 / 2 div (g)\\ Lq(Ma) < C g \\g\\ Lqmn , 2-e^q^2*+e. (4.97) 

Here, e > is a small number, not necessarily the same in each occurrence in 
(4.95)-(4.97), but which depends only upon n, a\ 1 and 02- We now claim that for 
some e > 0, we have 

\\(I + t 2 L a y^ 2 t 2 R(I + t 2 L a )-^ 2 \\ B{Lr{Rn}} «C C r e , 2-e<r<2 + £ (4.98) 
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(where of course C r depends also upon n, 01, and 02). Let us take this claim for 
granted momentarily. 

To prove (4.10), we first take p = q = 2 in (4.95), and r = 2 in (4.98), to obtain 
that for Eq < C2, the Neumann series in (4.94) converges in I? . We therefore obtain 
the case p = 2 of (4.10). The remaining cases of (4.10) then follow by interpolation 
from the endpoint case p — 2* — e, where we choose e so that (4.95) holds, and 
such that for this p, the exponent r := p* = pn/(n — p) < 2 lies within the range 
specified in (4.98). Let us then momentarily fix this choice of p. We then obtain 
that the Neumann series converges as a mapping from L p —> L 2 , by first applying 
(4.95) with q = p*, then (4.98) with r = p* , and with e < C p *, and then (4.95) 
with q = 2 and with p replaced by p* . This concludes the proof of (4.10) when 
n ^ 3. 

Similarly, we may obtain (4.12), and (4.13), in the former case by invoking (4.95) 
with p = q = 2, (4.98) with r = 2 and e Q < C 2 , and (4.96) with p = 2. To prove 
(4.13), we apply (4.97) with q = 2 - e, then (4.98) with r = q = 2 - e, and with 
Sq < C q , and then (4.95) with p replaced by the current q = 2 — e, and with q 
replaced by 2. 

To complete the proofs of (4.10), (4.12), and (4.13), with n 3, it remains to 
verify the claimed bound (4.98). To this end, let g £ L r (R"), and write 



\(I + t 2 L a )- 1 l 2 t 2 R{I + t 2 L a y 1 / 2 g\ 

+ t 2 L a )-^ 2 t 2 B x ■ v (/ + t 2 L a y 



1/2 



9 



L r (TR™) 



+ II (/ + t 2 ia)" 1/2 t 2 div (B 2 (I + t 2 L a )-V 2 g) || 

+ + t 2 L a )^ 2 t 2 V (I + t 2 L a )-^ 2 g\\ Lr{Rn) 

=: I + 11 + III. (4.99) 

By (4.95) with q = r, and p = r* := rn/(n + r), followed by Holder's inequality 
and (4.3), and then (4.96) with p = r, we obtain that 

I<C r eo\\g\\ L r<r), (4-100) 

as desired. Similarly, we may handle terms II and 177 by use of (4.95) and (4.97), 
with appropriate choices of exponents, along with Holder's inequality and (4.3). 
We omit the details, which are by now familiar. 

Let us now verify (4.11), (4.12), and the case p = 2 of (4.10), when n = 2. In 
this case, we use the resolvent expansion 

00 

(I + fMaY 1 =(I + t'La)- 1 J2[- f2 B- V(J + t'La)- 1 ] " , (4.101) 

fc=0 

where now div (£>) = 0. We note that {I + t 2 L a )~ 1 satisfies the following estimates 
in M 2 : 

\\{I + tL a ) 118(^1^2)^2^2)) + + t Lg) || B ( J3 -l( R 2) ) J,2( R 2)2) ^ C t , 

(4.102) 

and 

\\(I + tL a ) || B ( I/ 2( R 2) ii 2( R 2)) + ||*V(7 + t L a ) || B ( L 2( K 2) i 2( R 2)2) 

+ ||i 2 V(7 + i 2 J L a )- 1 div|| B(i2(R2)2jJ . 2(K2)2) sCC, (4.103) 
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Indeed, (4.103) is completely standard, and estimate (4.102) then follows by two- 
dimensional Sobolev embedding. By (4.7), (4.102), and a simple induction argu- 
ment, we have that 



(/ + fLa)- 1 [-t 2 B- V(7 + t 2 L a )~ x ] k f 



< r x C h e\ \\f\\x^M?) , (4.104) 



whence (4.11) follows immediately from (4.101), for Eq small enough. Similarly, for 
k ^ 1, using (4.7) and both (4.102) and (4.103), we obtain 

(7 + ^La)- 1 [-t 2 B- V(J + fLa)' 1 ] "J 



tV(J + t 2 L a )- x [-t 2 B- V(7 + t 2 L a )- 1 } k f 
< C^eg-^tT?- V(I + t 2 L a )- 1 f\\ H < C 



o 



L 2 (R 2 ) 2 
ll/IU 2 (I 



2 )- 



(4.105) 



By (4.101), we then obtain (4.12), and the case p = 2 of (4.10). The remaining 
cases 1 < p < 2, of (4.10) in two dimensions now follow by interpolation. 

Next, we turn to the proof of (4.15). The proof follows that of [13, Lemma 2.1], 
and is essentially the same in all dimensions n 2. To be precise, we give the proof 
in the case n 3, and then discuss, at the end of the argument, the modifications 
needed to treat the two-dimensional case. 

Remark 4.8. We note that M a has the property that its adjoint is of precisely 
the same class, thus, any estimates that we establish for M a can be immediately 
dualized. This property holds even in two dimensions, since in that case we have 
supposed that B is divergence free. 

We fix two cubes Qi,Q2, set d := dist(Qi, Q2), let r := mindiam((5i, Q2), 
R := maxdiam(<2i, Q2) and suppose that < t ^ 4r, and that 7? ^ 100(7. Let 
T] G C^°(IR n ), supported on a d/2 neighborhood of Q2. Thus, in particular, r\ = 
on Q\. We now set u := (7 + i 2 il7 a ) -1 /, where supp (/) C Q\. We claim that 



/^ |t v„ lV + /^ MV <^/^HW 

We verify the claim as follows. Observe that by definition 



d n x{I + t 2 M a )uur] 2 = / d n xfurf = 0, 



(4.106) 



(4.107) 



since / and 77 have disjoint supports. Rewriting the left-hand side of this last 
identity (using the definition of M a ), integrating by parts, and shuffling terms, we 
obtain 



t- 



d n xaVu-Vuri 2 + I d n x\u\ 2 r) 2 



d n xB\-VuuT] 2 + t 2 I d n xuB 2 -VuT] 2 - t 2 I d n xV\u\ 2 r] 2 



2t 2 J d n xa\7u-V7]ur) + 2t 2 J d n x B 2 \u\ 2 X7r]r] 



=: 7 + 77 + 777 + IV + V. 
By ellipticity, this implies that 



J d n x\tVu\ 2 ?f + J d n x\u\ 2 r, 2 < \I + II + III + IV 



V\ 



(4.108) 
(4.109) 
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Using Holder's inequality (equivalently, (4.27) and the dual version of (4.28)), and 
then Cauchy's inequality, we have 

2/2*\ 



\I + II + III\ < e t 2 ^ J d n x |Vu| V + (^j d n 



X \UT}\ 2 



<e t z (J d"x|V«|V + J d n x\V(urj)\ 

< s t 2 f^j d n x |Vu|V + J d n x \uVi]\ 2 ^j , (4.110) 
where in the next-to-last step we have used Sobolev embedding. Similarly, 

2/2*\ 



\V\ <£oi 2 I J d"a;|u| 2 |V?7| 2 + (j d n x\w 1 \ r 



< e t 2 I J d n x\Vu\ 2 if + J d n x\uy^\ 2 j . (4.111) 

For eo small enough, we may hide the integral involving Vw (which, by virtue of 
our previously established global estimate (4.12), is finite) on the left-hand side of 
(4.109). 

By Cauchy's inequality with e's, we have that 

\IV\<et 2 J d n x\Vu\ 2 V 2 + h 2 J ' cTx\u\ 2 \Vti\ 2 , (4.112) 

where e > is at our disposal. Taking s small enough, depending only upon 
ellipticity, we may, as before, hide the integral involving Vu on the left-hand side 
of (4.109). Collecting our various estimates, we obtain (4.106). 

Let us now deduce (4.15) as a consequence of (4.106). Recall that rj is supported 
on a d/2 neighborhood of Q2, which we denote by Q, and observe that diamQ ~ d 
(since, in particular, d > diam(Q2))- Following essentially verbatim the proof of 
[13, Lemma 2.1], we apply (4.106), with 77 of the form 

i] := e aip - 1, (4.113) 

with a := cod/t, where cq is a sufficiently small constant to be chosen momentarily, 
and where < ip € Cq°(Q), with ip = 1 on Q 2 , satisfying ||V^>||oo < V^- With this 
choice of rj, and with Co chosen sufficiently small, (4.106) yields in particular that 

f d n x \u\ 2 {e av - l) 2 s$ 1 J d n x \u\ 2 \e a *\ 2 . (4.114) 

Thus, by a simple argument involving the triangle inequality, we obtain in turn 
that 

P d n x\u\ 2 (e av ) 2 sC4 f d n x\u\ 2 <[ d n x\f\ 2 , (4.115) 



JI" JQi 

where in the last step we have used the case p — 2 of (4.10), and the fact that / is 
supported on Q\. Since ip = 1 on Q2, we then have that 

e 2a f d n x\u\ 2 < f d n x\f\ 2 , (4.116) 
JQi JQi 

whence the desired bound (4.15) follows for (I + t M a ) . 
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Next, let Q' 2 be a cube concentric with Q 2 , having side length £(Q' 2 ) ~ d, such 
that Q 2 cQ' 2 C Q, and such that dist(Q 2) M™\Q) « d « dist(Q 2 , R"\Q' 2 ). We note 
that the preceding argument may be applied with Q 2 i R place of Q 2 (indeed, we 
may construct ip enjoying the same properties as above, but with ip = 1 on Q' 2 ), so 
in particular (4.116) holds for Q' 2 . Thus, 

d n x\u\ 2 < e" cd /* [ d n x\f\ 2 , (4.117) 
Q'a JQi 

by definition of a. We now choose another cut-off ip e Cq°(Q' 2 ), such that ||V^||oo < 
1/d, with ^ (p ^ 1, and (p = 1 on Q 2 . We apply (4.106) with 77 = to obtain in 
particular that 



J d n x\tVu\ 2 ^ J d n x\tVu\ 2 (p) 2 <t 2 J d n x\u\ 2 \Vp\ 2 < (^j J d n x\u\ 2 

d n x\f\ 2 , (4.118) 



< p-cd/t 



where in the last inequality we have used (4.117) and the fact that t < d. We 
therefore have that (4.15) holds for tV(/ +t 2 M a ,) _1 , and thus by a duality argument 
(cf. Remark 4.8 above), it holds also for t(I + t 2 M a ) div(-). This concludes the 
proof of Lemma 4.4, when n ^ 3. 

The proof in the two-dimensional case is somewhat different: we shall use a 
resolvent expansion as in the proof of the global estimates treated above. This 
approach does not seem to yield the exponential decay that we obtained in higher 
dimensions, but nonetheless, as noted previously, polynomial decay of sufficiently 
high order suffices for our purposes. Let us first record some known estimates 
for the resolvent (J + i 2 L a ) _1 in two dimensions, which follow readily from two- 
dimensional estimates for the heat semigroup e~ sL proved in [18], and the formula 
(I + L)- 1 = J °° ds e~ s e~ sL . Let E, F C R" be disjoint sets, with 

d := dist(£,F) ^ t > 0. (4.119) 

We have 



\\(I + tL a ) || B ( £ l( £ ;) ii 2( f )) + ||*V(7 + t L a ) \\b(L 1 (E),L 2 (F)^) 

< Ct-^e-^l* < C d^it/d) 4 , 

and also 

|| (I + t 2 L a ) 1 \\ B{L 2 {E) L 2( F )) + + t 2 L a ) 1 || i3 ( L 2 (B) L 2 (F)n) < Ce 

< C{t/d)\ 



(4.120) 



(4.121) 



where as usual the various constants depend only upon a\ and a 2 . We now claim 
that there is a constant Co such that for / supported on E, 



H(k) := {I + t 2 L a )- 1 [-t 2 B-V{I + t 2 L a )- 1 ] k f 

^C fc ^'(t/d) 4 ||/|| L 2 (B) , 



L2 ( p ) (4.122) 
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and 



v(k) := tV(7 + i 2 i a ) _1 [ - t 2 S • V(7 + t 2 ^) -1 ] V 
< C%e k (t/d) 4 \\fh HE ). 



L 2 (_F)" 



(4.123) 



Given the claim, and choosing £o < 1/Co, we obtain (4.16) immediately from 
(4.101) and Remark 4.8. 

We now establish this claim by induction. Observe that the case k = is simply 
(4.121). Suppose now that (4.122) and (4.123) hold for some k ^ 0, for every set 
E and F as above, with d t, and we wish to verify the corresponding estimates 
for k + 1. Let 7i, F-i denote, respectively, d/3 and 2d/3 neighborhoods of F, and 
choose smooth a smooth cut-off ip such that 



= 1 on 7\ , supp (</?) C F 2 , 



(4.124) 



With ||V^|Uo 



*)' 



< 



1/d. We may suppose that t ^ d/3, since the case d ^ t > d/3 



follows from the global bound (4.105). We then have, for instance, that 
H(k + 1) = (I + t t L a )- 1 (-t l B + t*L a )~ l ) 



x [-t 2 B-V{I + t 2 L a )- l ] k f 
(I + fLa)- 1 {-t 2 B- V((p(I + fLa)- 1 )) 



2r \-l> 



L?(F) 



x [-t 2 B-V(I + t 2 L a )- 1 ] k f 

- I (J + i 2 ^)-^ - t 2 S • V((l - if){I + PL.)- 1 )) 

x [-t 2 B-V{I + t 2 L a )- 1 ] k f 
I + 11. 



L 2 (F) 



(4.125) 



By (4.120), since 1 - ip is supported on £1 := R n \7i, with dist(£i,F) = d/3, we 
have that 



77 < C(t/d) 5 iB • V((l - <p)(I + 1 2 ^)- 1 ) [-t 2 B- V(I + PLa)- 1 ] k .f 

< <7e (f/d) 5 |*V ((1 - p)(J + t 2 ^)- 1 ) [ - t 2 B ■ V(7 + t 2 La) _1 ] */ 

< Ce (t/d) 6 | i 1 + **L a )- 1 [-t 2 B- V(7 + t 2 L a )- 1 ] k f" 
+ Ce (t/d) 5 \\tV(I + t 2 L a )- 1 [-t 2 B- V(7 + t 2 L a )~ l ] K f 

<C k+1 e k +1 (t/df\\f\\ LHE) , 



L 1 (R™\F 1 ) 
L 2 (R™)™ 



L 2 (R") 



(4.126) 



where in the second inequality we have used Holder's inequality and (4.3), in the 
third inequality that |[V^||l°°(R'»)'> < and in the last inequality we have used 
the global bound (4.105). 

Turning to term 7, we have by (4.102) and (4.7), since supp (ip) C F2, with 

I|Vp||l-o 



< 



7 < Cec 



1 /d, that 

tV(cp(I + t 2 ^)- 1 ) [-t 2 B- V(7 + t 2 7 a )-!] fe / 



i 2 (F 2 )-- 
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k 



< Ce *V(J + 1 2 ^)- 1 [-t 2 B- V(J + ^La)" 1 ] 7 

+ Ceo (t/d) J (7 + 1 2 ^)- 1 [ - t 2 B • V(J + 1 2 ^)" 1 ] */ 
^2CC fc eg +1 (3t/d) 4 ||/|| i2(£) , 



(4.127) 



by the induction hypothesis, since dist(i r 2, E) d/i. Combining our estimates for 
terms I and II, we obtain the case k + 1 of (4.122), for a suitable choice of Co- 
By a very similar argument, one may establish the case k + 1 of (4.123). We omit 
further details. □ 



Appendix A. Additional Facts on Kato's Resolvent Approach 

In this appendix we further investigate T as defined in (2.4) and (2.5) and also 
make the connection between T and the form sum of To an d B* A under stronger 
hypotheses on T , A, and B compared to those in Hypothesis 2.1. 

For relevant background literature in this connection we refer, for instance, to 
[34, Ch. IV], [42], and [67] (see also [39], [84], [85, Ch. II] in the context of self- 
adjoint operators). 

We start by briefly complementing (2.4): 

Lemma A.l. Assume Hypothesis 2.1. If T is defined by (2.5), then T* is given 
in terms of Tq , A, and B by 

(T* - zlnY 1 = (T * - zlnY 1 

- (T*-zI H )-^A*[l K - K(z)] 1 B(T^ - zI H y\ (A.l) 

ze{(e P (T*)\iep(K(0)}, 

where 

K(z) = —B(Tq* — zI- H )~ 1 A* 1 z G p(T*). (A.2) 

Therefore, T* is defined via Kato's resolvent equation by replacing To by Tq and by 
interchanging the roles of A and B. 

Proof. Taking adjoints on both sides of (2.6), and making use of the identities 

[A(To - zlu)- 1 ]* = [(To* -zI H )- l A*)*}* = (T Q * -zlnY'A* (A3) 

[{To-zI H )-^B*\* = [(To - zIur'B*]* = B(TZ-zI n )~\ (A.4) 
one obtains the following representation of the resolvent of T Q *: 
(T^-zIu)- 1 

= (TZ-zIu)- 1 - [A(T - zlu)- 1 } * {[lie - K^)]- 1 )* [(T - zI n )-^B*} * 

= (T* - zl n )- x - (T* -zI n )-iA*(I K - K(z)*)~ 1 B(T - zI H )-\ (A.5) 

z e {( e p(T a )\i e P (K (())}. 

Finally, a short computation with adjoints and closures shows 

K(z)*=K(z), ze P (T ). (A.6) 
Thus, (A.5) and (A.6) combine to yield (A.l). □ 
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Before stating the principal result of this appendix, Theorem A. 3 below, we 
need to recall some preparations. The material needed can be found in [34, Sect. 
IV. 2], [37, Sects. 1.2-1.5], and [68, Sect. VI. 2. 6], and we here follow the summary 
presented in [43, Sect. 2], except, the latter focused exclusively on the self-adjoint 
situation, whereas here we extend it to the non-self-adjoint (sectorial) context: Let 
V be a reflexive Banach space continuously and densely embedded into %. Then 
also "H embeds continuously and densely into V*. That is, 

V^n^V*. (A.7) 

Here the continuous embedding H ^-y V* is accomplished via the identification 

H 3«4 (-,v)n eV*, (A.8) 

and we use the convention in this manuscript that if X denotes a Banach space, 
X* denotes the adjoint space of continuous conjugate linear functionals on X, also 
known as the conjugate dual of X. 

In particular, if the sesquilinear form 

v( v )v.:VxV*4C (A.9) 

denotes the duality pairing between V and V*, then 

v(u,v) v * = (u,v) n , ueV,veH^V*, (A.10) 

that is, the V, V* pairing v( ' > ' )v* is compatible with the scalar product ( • , • 
vaH. 

Let W G S(V,V*). Since V is reflexive, (V*)* = V, one has 

W:V^V\ W*:V^V* (A.ll) 

and 



v {u,Wv) v , = v .(W*u,v) {v * r = v *(W*u,v) v = v (v,W*u)v*. (A.12) 
Self-adjointness of W is then defined by W — W* , that is, 



v{u, Wv)v = v* (Wu, v) v = v (v, Wu) v * , u, v G V, (A. 13) 
nonnegativity of W is defined by 

v (u,Wu) v * ^0, u e V, (A.14) 
and boundedness from below of W by c\y G R is defined by 

v (u,Wv)v >c w \\u\\^, ueV. (A.15) 

(By (A. 10), this is equivalent to v(u,Wu)v* > cw v( u , u )v* ■ u G V.) 

Next, let the sesquilinear form o( • , ■ ) : V x V — > C (antilinear in the first and 
linear in the second argument) be V-bounded, that is, there exists a c a > such 
that 

|a(u,v)|<c a ||u||v|M|v, u,v€V. (A.16) 

Then A defined by 

A: { . . (A.17) 
I v t-y Av = a( ■ , v), 

satisfies 

A eB(V,V*) and v (u,Av) =a(u,v), u,veV. (A.18) 
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Assuming further that o( • , • ) is symmetric, that is, 

a(u, v) = a(v, u), a.ueV, (A. 19) 
and that a is V -coercive, that is, there exists a constant Co > such that 

o(u,u)>Cb||«||v > ueV, (A.20) 

respectively, then, 

A: V — > V* is bounded, self-adjoint, and boundedly invertible. (A. 21) 
Moreover, denoting by A the part of A in H defined by 

dom(A) = {u e V \ Au e %} c H, A = A\ Aom , A) ■ dom(A) -> H, (A.22) 
then A is a (possibly unbounded) self- adjoint operator in H satisfying 

A > C I n , (A.23) 

dom {A 1 ' 2 ) = V. (A.24) 

In particular, 

A- 1 e B(U). (A.25) 

The facts (A.7)-(A.25) are a consequence of the Lax-Milgram theorem and the 2nd 
representation theorem for symmetric sesquilinear forms. Details can be found, for 
instance, in [32, Sects. VI.3, VII.l], [34, Ch. IV], and [71]. 

Next, we extend the scope of this discussion and consider a sectorial form 
b( - , ■): V x V — > C defined by the existence of c b E K and 9 b E [0, 7r/2) such 
that 

{b(u, u) E C | u E V, \\u\\ H = 1} C {C G C | | arg(C - c 6 )| < b }. (A.26) 
Denoting 

bfn(u, v) = Re(b(it, v)), u, v E V, (A. 27) 

this yields 

b m (u,u) ^ c b \\u\\ 2 H , ueV. (A.28) 

Introducing the scalar product ( ■ , • )v : V x V — > C (and the associated norm 
II -Ik) by 

(u,v) Vb = b<x{u,v) + (1 - c b )(u,v) H , u,veV, (A.29) 

turns V into a pre-Hilbert space (V; ( • , • )y ), which we denote by V b . The form 
b is called closed in TL if V b is actually complete, and hence a Hilbert space. The 
form b is called closable in % if it has a closed extension in H. The estimates, 

\bm,(u,v) + (l-c b )(u,v) H \^ HvJMIv,. u,veV, (A.30) 

and 

\bn{u,u) + {l-c b )\\u\\l\ = \\u\\\, ueV, (A.31) 

show that if b is closed in then the form b^( ■ , • ) + (l — c b )( ■ , •)% is a symmetric, 
V-bounded, and V-coercive sesquilinear form and hence the analog of (A.17)-(A.25) 
applies. Similarly, one has the bound 

|b(u,i>) + (l -<*)(«, «) w | < [l + tan(ff 6 )]||«||vj|«||v b , u,veV. (A.32) 
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Hence, if b is closed, the 1st representation theorem for forms (cf. [34, Theorem 
IV. 2. 4], [68, Theorem VI. 2.1]) yields the existence of a linear map 

B ■ < b ^' (A 33) 

C "' \v^B Cb v = b(-,v) + (l-c b )(-,v) H> 

with 

B Cb eB{V b ,V;) and Vb (u,B Cb v) v , = b(u, v) + (l-c b )(u, v) H , u,v eV. (A.34) 

Introducing the linear map 

B = B Ci + (c& — 1)1: V b -> V£, (A.35) 

where / : Vb V£ denotes the continuous inclusion (embedding) map of V b into 
V ( *, one obtains an m-sectorial operator B in % by restricting _B to 

dom(B) = {u G V | Bu G U) C "H, B = £| dom(B) : dom(S) -> H, (A.36) 
satisfying the following properties: 

b(u, v) = Vb (u,Bv) v ,, u,veV, (A. 37) 

b(u, v) = (u, Bv) H , u G V, v E dom(B), (A.38) 
dom(i?) = {v G V | there exists an /„ G such that 

E>0, v) = (to, for all w G V}, (A.39) 
Bii = / u , u G dom(_B), 

dom(i?) is dense in T-L and in V b . (A. 40) 

Properties (A.39) and (A. 40) uniquely determine B; the operator B is called the 
operator associated with the form b; conversely, the form b is called the form asso- 
ciated with B. 

Finally, we also recall the following perturbation theoretic fact: Suppose that a 
is a sectorial and closed form with domain VCW, and let b be a sesquilinear form 
bounded with respect to a with bound less than one, that is, dom(b) 3 V, and that 
there exist ^ a < 1 and j3 > such that 

\b{u,u)\ ^a\a{u,u)\ + P\\u\\^, ueV. (A.41) 

Then 

I (u, v) (a + b)(u, v) = a(u, v) + b{u, v) 

defines a sectorial and closed form (cf., e.g., [68, Sect. VI. 1.6]). In the special case 
where a can be chosen arbitrarily small, the form b is called infinitesimally form 
bounded with respect to a. In the case where T a is the m-sectorial operator uniquely 
associated with the form a, and also b is sectorial with uniquely associated m- 
sectorial operator T bl the m-sectorial operator uniquely associated with the form 
sum o + b in (A. 42) will be denoted by T a+b as well as by T a + q T b , that is, we 
agree to use the notation T a+b = T a + q T b in this situation. 

Given these preliminaries, we now turn to the connection between Kato's resol- 
vent equation (2.4) defining T and the form sum of Tq and an extension of B* A 
under somewhat stronger hypotheses on Tq, A, and B than those in Hypothesis 
2.1. 

We start with the following basic assumptions: 
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Hypothesis A. 2. (i) Suppose To is m-sectorial and 

dom (T 1/2 ) = dom ((T*) 1 ' 2 ). (A.43) 

(ii) Assume that A : dom(A) — > IC, dom(A) CM, is a closed, linear operator from 
TL to IC, and B : dom(_B) — » IC, dom(i?) CM, is a closed, linear operator from T-L 
to IC such that 

dom(A) D dom (T 1/2 ) , dom(B) D dom (T 1/2 ) , (A.44) 
and that for some ^ a < 1, b 0, 



IIA/H 2 ; ^ aRe[q To (/,/)] +6]]/]|^ / € dom(T 1/2 ), 
\\Bf\\l < aRe[q To (/,/)] +&||/|&, / £ dom(T 1/2 ). 



(A.45) 



(Hi) Suppose that 



where 



lim \\K(-E)\\ m) =0, (A.46) 



X (z) = -A(T - zI n )-^B* 



= -[A{T a -zI H )- 1/2 } [B(T* -zln)' 1 ' 2 ]* 6 B(K), z e P (T ). 



(A.47) 



In the following we denote by q T the sectorial form associated with To and note 
that according to Remark 2.16, assumption (A.43) yields that 

dom (T 1/2 ) = dom(q To ) = dom ((T*) 1 / 2 ). (A.48) 

In addition, introducing the form 

q w (f, g) = (Bf, Ag) K , f, g G dom (T 1/2 ) (A.49) 

(formally corresponding to the operator B*A), it is clear from assumption (A.45) 
that q w is relatively bounded with respect to q Tg with bound strictly less than one 
(cf. (A. 41)). Thus, we may introduce the sectorial form 



qr 0+q w(/»ff) = qT o (/>0) + 1w(f,9), 
f,g G dom(q To+qH/ ) = dom(q T J = dom(T i/2 ), 



1/3 (A.50) 



and hence denote the m-sectorial operator uniquely associated with q To+ ^ w by 
Tq+qW in the following. 

The principal result of this appendix, establishing equality between T defined 
according to Kato's method (2.4), (2.5) on one hand, and the form sum To + q W 
on the other, a result of interest in its own right, is proved next: 

Theorem A. 3. Assume Hypothesis A. 2. Then T defined as in (2.4) and (2.5) 
coincides with the form sum of Tq + q W in T-L, 

T = T + q W. (AM) 

Proof. The strategy of proof (closely following the one in [85, Theorem 11.34] in 
connection with Rollnik potentials on R 3 ) is to show that the resolvents of T and 
To + q W both coincide with the right-hand side of (2.4). 
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For this purpose we introduce V T and V£ in analogy to V b and V£ with the form 
b replaced by q To in (A.26)-(A.32), and similarly the extension To : V Tq — > V^ o of 
To : dom(To) — > H as in (A. 35). Moreover, since for some fixed C > 0, 

< \\Bf\\ K \\Ag\\ K <C\\f\\ VTo \\g\\ VTo , /, g G dom (T Q 1/2 ) , (A.52) 

there is a natural map 

WeB(V To ,V£ ), 

, 53) 

(Wff) (/) = v rn </, W 5 > v , = (-B/, A fl ) K) /, g G V To , 

extending 5*74, such that 

T ~M/ = 7b + G B(V To , V£ ) , 

{{n+^W)g) (f) = (f g) (/) + (W<?) (/) = v </, Tbff) v . + V r „ </, W ? 9> v . 
= qT (/. 5) + A fl ) K , /, .g G V To . (A.54) 
In addition, we note that B* : dom(B*) — » H, dom(_B*) C /C, extends to 

B* eB(JC,V^ ), 

(B*k)(f) = VTo (f,B*k) v ^ = ((To + c/«) 1/2 /, [B(T + cl n )- 1/2 ]*k) n , (A.55) 

/g V To , ke!C. 

This implies 

W = B*A, (A.56) 

and also 



K(z) = -A(To-zI u )- 1 B* = -A(To-zl) B* G B(IC), z G p(T ), (A.57) 
where A : dom(A) — > /C, dom(A) C H and 

(Tb - zTy 1 G B(V£ , V To ), z G p(T ). (A.58) 
Next, we introduce the operator 

R(z) = (To - zlnY 1 -(To- zl)~ l B*[I K - K(z)]- 1 A(T a - zI H T\ (A gg . 

z G {C G p(T ) | 1 G p(K(C))}, 

and hence conclude that 

R(-E) = (T + EIu)- 1 (A.60) 

- (r + T/)" 1 B*^(-l) fe [A(To + T/)" 1 B*J A(To + EI H )-\ 

k=0 

for £7 > sufficiently large such that 

\\K(-E)\\ B w < 1. (A.61) 

We note that the infinite sum in (A. 60) is convergent in Z?(/C) and hence the mapping 
properties in (A.53)-(A.58) yield R(-E) G B(U) for E > sufficiently large. One 
then computes, taking into account telescoping of series, that 



(To + EIuyHTo+^W + EI] 



(To + EIu)- 1 
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N 



-(T + EI) B*J2(-l) k [MTo 

k=0 

(T + EI H )- 1 (lb+B'*A + EI 

All 



N 

(T +Eiy^B*^2(-l) 

k=0 

= {t + ei h )- 1 

+ (-l) N+1 {T Q + EI 

iVfoo 



EI) 1 B* 



(T + Ei H y 



EI) B* 



A{T + EI n ) 



A(T +EI n ) 



B* 



A To + EI 



B* 



N+l 



A(T + EIn)' 



(T + EIn)- 1 in £(^)-norm. 



Thus, 



(T + Eln)- 1 {T + B*A + EI)R(-E) = (Tq+EIu)- 1 
= (T + EIn)' 1 (To + B*A + El) (T + B*A + El) ' 



implying 



R(-E) = (T + B*A + EI ) for E > sufficiently large, 



(A.62) 



(A.63) 



(A.64) 



since ker (T + El) = ker (T + B* A + El) = {0} for for E > sufficiently large. 
Analytic continuation with respect to E then proves 

-l 



R{ z ) = [To + B* A - zl 



z G {( 6 p(T ) 1 1 G p(K(C))}. 



(A.65) 



A comparison of the right-hand sides of (2.4) and (A. 59) then yields R(z) = R(z) 
and hence (A. 51). □ 

Acknowledgments. We arc indebted to Alan Mcintosh for helpful discussions. 
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